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Previous work on electroweak radiative corrections to high energy scattering using soft-collinear 
effective theory (SCET) has been extended to include external transverse and longitudinal gauge 
bosons and Higgs bosons. This allows one to compute radiative corrections to all parton-level hard 
scattering amplitudes in the standard model to NLL order, including QCD and electroweak radia- 
tive corrections, mass effects, and Higgs exchange corrections, if the high-scale matching, which is 
suppressed by two orders in the log counting, and contains no large logs, is known. The factorization 
structure of the effective theory places strong constraints on the form of gauge theory amplitudes 
at high energy for massless and massive gauge theories, which are discussed in detail in the paper. 
The radiative corrections can be written as the sum of process-independent one-particle collinear 
functions, and a universal soft function. We give plots for the radiative corrections to qq — > WtWt, 
ZtZt, WlWl, and ZlH, and gg — > WtWt to illustrate our results. The purely electroweak cor- 
rections are large, ranging from 12% at 500 GeV to 37% at 2 TeV for transverse W pair production, 
and increasing rapidly with energy. The estimated theoretical uncertainty to the partonic (hard) 
cross-section in most cases is below one percent, smaller than uncertainties in the parton distribu- 
tion functions (PDFs). We discuss the relation between SCET and other factorization methods, 
and derive the Magnea-Sterman equations for the Sudakov form factor using SCET, for massless 
and massive gauge theories, and for light and heavy external particles. 



I. INTRODUCTION 

Soft collinear effective theory (SCET) 0, 0, S, H is a 
field theory which describes the interactions of energetic 
particles which produce final states with small invariant 
mass. SCET was developed for QCD, but has recently 
been extended to spontaneously broken gauge theories 
with massive gauge bosons [E IE 0- This allows one 
to compute electroweak corrections to Standard Model 
processes at high energies, and to sum electroweak Su- 
dakov logarithms [E0|- One can compute QCD and elec- 
troweak corrections to parton-level scattering amplitudes 
with 1% precision. To achieve this accuracy, one needs 
to compute two-loop QCD and one-loop electroweak ra- 
diative corrections, and sum the large logarithms using 
renormalization group methods in SCET. An extensive 
analysis with plots of numerical results for 2 —> 2 parton 
processes such as dijet production, Drell-Yan, and top- 
quark production was given in Ref. Q. The results in- 
cluded fermion mass effects (e.g. due to the top quark), 
7 — Z mixing (which leads to Myy ^ Mz) as well as 
Higgs exchange corrections. The results are extended in 
this paper to include Higgs and gauge bosons. The purely 
electroweak corrections are substantial — 20% for W pair 
production and 10% for ZH production at 1 TeV rising 
to 50% and 25% respectively, at 4 TeV, and are much 
larger than two-loop QCD corrections. 

The Sudakov double-logarithms arise as a result of soft 
and collinear singularities in radiation, and are a fea- 
ture of exclusive amplitudes such as form-factors. They 
are absent in totally inclusive processes such as the total 
e + e~ hadronic cross-section. There is some cancellation 
of the QCD Sudakov corrections in inclusive processes 
such as jet rates, because the infrared singularities can- 



cel between real and virtual graphs. Jet cross-sections 
still have a residual double-logarithmic dependence on in- 
frared parameters in the jet definition such as y cu t, and 
can have large radiative corrections. Unlike QCD, the 
electroweak Sudakov corrections do not cancel, even for 
a totally inclusive cross-section, because the initial state 
particles are not electroweak singlets [H, IE [l(| ■ 

Radiative corrections to high energy processes have 
been obtained previously usin g fix ed order computations 
by manygroupsjiji M. ELM SE, Ml EE El El El, 
EI HE Si M, Hi TIElllTll ■ The compu- 

tations use infrared evolution equations [Tlj |. based on 
an analysis of the infrared structure of the perturbation 
theory amplitude and a factorization theorem for the Su- 
dakov form factor [3(| ■ We have checked our EFT results 
against these computations. 

The computation of electroweak radiative corrections 
is greatly simplified by the use of effective field the- 
ory methods. The computation splits into a high-scale 
matching condition, and the EFT computation. The 
matching can be computed neglecting all low energy 
scales such as fermion masses, gauge boson masses, and 
spontaneous symmetry breaking. Secondly, the effec- 
tive theory results depend only on the color structure 
of the amplitude, and are independent of the momen- 
tum structure. This is because the large momenta (i.e. 
the kinematic variables s,t,u) are labels in the effective 
theory, and not dynamical variables. Thus an amplitude 
A(s, t, u) in the effective theory can be treated as a num- 
ber, with a different value depending on the kinematics. 
In contrast, in a computation without the use of effective 
theory methods, A is a vertex in a Feynman graph, and 
the arguments of A get integrated over in loop graphs. In 
gauge boson production, there are 10 possible kinematic 
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structures [3l[, and they all evolve the same way in the 
effective theory. 

In this paper, we explain in more detail the factoriza- 
tion structure of the high energy amplitudes which fol- 
lows from the effective theory jEJ, 0, 0, [HJ, and extend 
our results to include external gauge bosons and Higgs 
scalars, so that we have the results for all standard model 
particles. The computations for longitudinally polarized 
gauge bosons use the Goldstone boson equivalence the- 
orem HI M, H H3, HI S [!}. The factorization 
results hold both for the anomalous dimension, as well as 
for the low-energy matrix elements in the effective theory, 
i.e. for both the infinite and finite parts of the amplitude. 

The factorization structure is rederived using the A- 
regulator introduced in Ref. [4lJ] , because this gives the 
answer in a form which can be most easily compared with 
previous results using other methods IH, 0, EH ■ We 
show that the hard scattering amplitude can be writ- 
ten as a sum of process-independent onc-particle collincar 
functions which depend on external particle energies, and 
a universal soft function which depends only on the exter- 
nal particle directions. The amplitude for a hard scatter- 
ing process with r external particles is given by combin- 
ing the collinear functions for each external particle with 
the universal soft function. The collinear function for 
all standard model particles, and the universal soft func- 
tion can be computed using the results given in this pa- 
per. The explicit expressions for the standard model are 
rather lengthy, because of the different possible fermion 
quantum numbers, and because of custodial SU (2) vi- 
olation which leads to 7 — Z mixing, Mw 7^ Mz, and 
to different radiative corrections for the W and Z. For 
this reason, in this paper, we give the results for a broken 
SU(2) gauge theory, with 2np massless fermion doublets. 
The standard model results can be obtained straightfor- 
wardly from the SU (2) results, and will be given explic- 
itly in a companion paper [46[ . All numerical plots shown 
here arc in the full standard model, including Higgs cor- 
rections due to the top-quark Yukawa coupling. 

The high energy behavior of scattering amplitudes for 
massive gauge theories is the same as that for mass- 
less gauge theories, since symmetry breaking is a soft 
effect and does not influence the high energy behavior 
of the amplitude. Thus our results are related to recent 
work on the infrared structure of gauge theory ampli- 
tudes 0, S3, IH, [H, M, HI HI, which we discuss in 
Sec. |Xj Of particular interest is whether the SCET am- 
plitudes obey two properties, Casimir scaling and the 
sum-on-pairs form. 

We compute the Sudakov form factor in a massive 
gauge theory, following the analysis of Magnea and Ster- 
man [12]. This clarifies the connection between SCET re- 
sults and the analysis using factorization methods. While 
SCET and factorization methods [53| use the same termi- 
nology of soft, collinear and hard amplitudes, the differ- 
ent contributions are not identical in the two approaches, 
even though their sum is the same. In our computation, 
the infrared sector of the gauge theory is under complete 



perturbative control, since the gauge bosons are massive, 
and the gauge theory is spontaneously broken. Thus we 
can compute the low-energy matrix elements, which are 
finite and free of infrared divergences. This provides ad- 
ditional insight into the structure of gauge theory ampli- 
tudes. In the usual QCD analysis, the low energy matrix 
elements are infrared divergent, and infinite in perturba- 
tion theory. 

Using the results in this paper, one can compute ra- 
diative corrections in the effective theory to all partonic 
hard scattering processes, and thus to any hard scattering 
process at the LHC with an arbitrary number of external 
particles, if the high scale matching (step (2) below) is 
known. The high-scale matching is suppressed by two 
powers of a large logarithm relative to the EFT correc- 
tions. In the examples we have studied, the high scale 
matching (including QCD) is much smaller (by a factor 
of four or more) than the electroweak EFT corrections, 
so one can profitably use our EFT results even for cases 
where the high scale matching is not known. The only 
exception is the cross-section for qq — > W^W^f , where 
the high-scale matching is almost 20%. The radiative 
corrections are also large for this case, and exceed 20% 
below 1 TcV. 

The computation of the complete hadronic cross- 
section can be done using the following steps, shown 
schematically in Fig. [TJ 

(1) If the amplitude is the matrix element of a full the- 
ory operator, such as the Sudakov form factor for the 
scattering of a gluon by the operator 0{(Xq) = G^G A ^ , 
scale the operator using the standard model anomalous 
dimension 70 to the scale \ih (high-scale matching). If 
the amplitude is a scattering amplitude computed using 
vertices in the standard model Lagrangian, this step is 
not necessary, since the Lagrangian is /i independent. 

(2) Compute the hard matching 1 coefficient C at a 
scale fih- This is a standard example of a matching com- 
putation to an effective field theory. C is computed from 
graphs in the full theory, setting all small scales such as 
the gauge boson and fermion masses to zero. The re- 
sulting graphs are single-scale graphs, and have ultravi- 
olet and infrared divergences, which are regulated using 
dimensional regularization, leading to 1/e terms. The 
ultraviolet divergences are cancelled by renormalization 
counterterms, and the infrared divergences are repro- 
duced by corresponding infrared divergences in the effec- 
tive theory. Thus C is given by keeping the finite part of 
the full theory diagrams, and throwing away all singular 
terms in 1/e. 2 The hard matching is independent of low 



1 pun intended 

2 The above procedure only works if dimensional regularization is 
used to regulate both the ultraviolet and infrared divergences. 
A check on the matching computation is that the 1/e infrared 
divergent terms must agree with the anomalous dimension in 
the effective theory. See, for example, Refs. [53 . [5jj[ for a more 
extensive discussion. 
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FIG. 1: Steps in the computation of a hadronic cross-section 
at high energy. The effective theory can be used for everything 
except the hard matching correction C, and the running 70 
in the full theory. 



scales such electroweak symmetry breaking and fermion 
masses, and so can be computed in the unbroken theory 
using gluons, W and B gauge bosons, without worrying 
about electroweak gauge boson mixing, even for processes 
such as Z production. All mass effects are included in the 
EFT computation in steps (4)-(6). 

All scattering amplitudes are independent of the choice 
of /i/j. However, in perturbation theory, there is residual 
fih dependence from higher order terms that are not in- 
cluded in the computation. \ih is chosen to be of order 
the hard scale in the scattering problem so that there are 
no parametrically large logarithms in C, and all large 
logarithms are summed by renormalization group evolu- 
tion. We will choose (if, — the partonic center of 
mass energy in the collision. 

The high-scale matching for processes involving a small 
number of external partons is known and is included in 
the examples considered in Ref. and in this paper. 
The high scale matching, in general, does not obey the 
factorization structure of the effective theory amplitude. 



This is the only piece of the computation which cannot 
be computed using the effective theory results. 

(3) Run the amplitude from the high scale fih to a 
low scale \x\ of order Mz using the SCET anomalous di- 
mension computed in SCETew- ( SCETew is SCET 
with dynamical SU(3) x SU(2) x U(l) gauge bosons). 
The SCET anomalous dimension is linear in log Q 2 / pi 2 to 
all orders (56j, and the renormalization group evolution 
sums the Sudakov double-logarithms. For Q ~ 1 TeV, 
L 2 = log 2 Q 2 /Mg ~ 21 so the Sudakov corrections can be 
very large, and lead to a suppression of the cross-section. 
It is well-known that the QCD Sudakov corrections can 
have a huge effect on the cross-sections. The electroweak 
Sudakov corrections are also significant, ranging from 
10%-50% depending on the process at energies of a few 
TeV, increasing rapidly with energy. The Sudakov double 
logarithms which are summed by the LL series are two 
powers of L more important than the high-scale matching 
in step (1). The anomalous dimension is independent of 
the electroweak scale, and like the high-scale matching, 
can be computed in the unbroken theory. 

(4) Compute the low-scale matching at fii by integrat- 
ing out the W, Z, Higgs and i-quark, onto a theory con- 
taining only QCD and electromagnetism,which we call 
SCET 7 . The new feature found in Ref. 0, 0| is that 
this low-scale matching in SCET contains a \ogQ 2 //j,f 
term of order L. This term is related to the cusp anoma- 
lous dimension. No higher powers of a log appear to all 
orders in perturbation theory. For most processes, the 
low-scale matching correction is purely electroweak, since 
the gluon graphs are continuous through the electroweak 
threshold. The low-scale matching D can be comparable 
to or slightly larger than the QCD contribution to the 
high scale matching C because of the log Q 2 /p 2 enhance- 
ment. There are low-scale QCD corrections for processes 
involving external i-quarks, and here the corrections are 
much larger, of order 35%. The low-scale matching is one 
power of L more important than the high-scale matching. 

(5) Run the amplitude to \if (the factorization scale) 
using the anomalous dimension in SCET 7 . ( SCET 7 is 
SCET with dynamical gluons and photons). This gives 
the parton level scattering amplitude at the scale /i/, i.e. 
the hard amplitude at /.t = [if, denoted by H in Fig. [1] 

(6) The last step depends on the experimental ob- 
servable being computed. To compute a hadronic cross- 
section, one would first square the parton scattering am- 
plitude and multiply by the partonic phase space to get 
the hard partonic cross-section er(/i/) oc \H(/if)\ 2 . An 
SCET amplitude with two incoming collinear fields and 
n outgoing collinear fields can be used to compute the 
cross-section for pp — > n jets. This is an inclusive cross- 
section at the level of partons, i.e. each jet can contain 
any number of collinear partons, but it is exclusive at the 
level of jets, i.e. one only includes exactly n jets, and no 
additional hard radiation. For such exclusive jet rates, 
one in general can have parton distribution functions 
(PDFs) describing the parton distributions inside the 
proton, as well as jet and soft functions (53[. These must 
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all be scaled to the common scale u / , and then convoluted 
with the partonic cross-section a (fif), to get the experi- 
mentally measurable hadronic cross-section. The /if de- 
pendence of &(/j,f) is cancelled by that in / <E> f <£> J <£> S. 
This provides non-trivial constraints on the anomalous 
dimension, since the convolution running of /, J and 
S must be cancelled by the multiplicative running of 
H [57], UH- For complete consistency, the parton dis- 
tribution functions need to be extracted including QED 
corrections. The QED running below Mz is very small, 
so one can simply drop this and match to PDFs obtained 
neglecting electromagnetic corrections. 

The final step (6) above is common to all hadronic com- 
putations, and has been studied extensively. We have 
nothing new to say here. We therefore stop our com- 
putation at the end of (5). In the numerical results we 
present, we will choose fif — Mz- 

Various software packages such as BlackHat [59| . 
Rocket [6(| and CutTools [6l[ are being developed 
to streamline the computation of QCD radiative cor- 
rections to experimentally relevant cross-sections at the 
LHC. These calculations typically include the QCD ef- 
fects discussed in (l)-(4) above, but not the electroweak 
radiative corrections. They also perform the difficult task 
of taking partonic cross-sections and computing infrared 
safe observables (such as jet rates) in proton scattering. 
The electroweak corrections are independent of the QCD 
corrections to one-loop, and the mixed QCD-electroweak 
terms from renormalization group running at two-loops 
are tiny. One can then incorporate our results in the 
existing methods by using the results in this paper with 
a s — > as a multiplicative correction factor to the scat- 
tering amplitudes. If the partonic amplitude is needed at 
a scale Uf > Mz, it can be obtained by using our results 
for a s — ► and \i$ = Mz, and then scaling back up to 
fif using the pure QCD expressions. 

The outline of the paper is as follows. In Sec. HH we 
give plots of a few representative processes to illustrate 
our results. We consider ulul — > clCl as an example of 
fcrmion scattering, ulul — > WtWt, ZtZt for transverse 
electroweak gauge boson production, Ulul — » WlWl for 
longitudinal gauge boson production, ulul — > HZl for 
associated Higgs boson production, and gg — ► WtWt- 
The kinematics and notation is given in Sec. IIII| and the 
structure of the logarithmic series in Sec. IIVI Sec. IVII 
discusses the EFT result that we need for our computa- 
tions, which have been checked by explicit computation. 
There follows a long discussion in Sec. IVIH -Scc. IXIII on 
the general form of the radiative corrections to all orders 
in perturbation theory, which can be skipped by readers 
not interested in technical details. The topics discussed 
are only relevant beyond three loops, and away from the 
N — > 00 limit. Sec. IVHI discusses the A-regulator, and 
Sec. IVIIll discusses the factorization structure of the am- 
plitude. The general form of the scattering amplitude 
is given in Sec. IIXI Casimir scaling, the sum-on-pairs 
form, and A"-factors are given in Sec.|X] The decomposi- 
tion into collinear and soft functions is given in Sec. IXI1 



and the extension to heavy quarks in Sec. IXIII The Gold- 
stone boson equivalence theorem, and the one-loop radia- 
tive corrections, are summarized in Sec. IXIIII Sec. IXIVI 
gives the results for the collinear and soft functions in an 
SU(2) gauge theory. The Magnea-Sterman result for the 
Sudakov form factor is given in Sec. IXVI and conclusions 
in Sec. IXVll 



II. PLOTS 

The numerical results for a few representative standard 
model computations are shown here. The details of how 
these plots were computed will be discussed later in this 
article. All the computations have been done using QCD 
plus electroweak one- loop matching at pb% and fxi, and 
renormalization group evolution using the SCET QCD 
plus electroweak non-cusp anomalous dimension to two 
loops, and the QCD plus electroweak cusp anomalous 
dimension to three loops. The electroweak three- loop 
cusp anomalous dimension is only known neglecting the 
Higgs scalar contribution. 

The results have a very small dependence on the Higgs 
mass. All the numerical results are for to# = 200 GeV. 
Changing tuh to 500 GeV changes the rates by less than 
one part in 10 4 . 

The full standard model is matched onto SCETew at 
fib, and the electroweak gauge bosons, Higgs, and top 
quark are integrated out at ui- We will choose Hh = 
y/I and Hi = Mz in computing scattering rates. It is 
convenient to integrate out the electroweak gauge bosons, 
Higgs, and top quark at a common scale, so that one does 
not have to treat an effective theory with electroweak 
gauge bosons but broken gauge symmetry. 

The fih,i dependence is very small in almost all the 
processes we have studied. Fig. [2] shows the variation in 
the ulul — * clCl cross-section as these scales are varied. 
Uh is varied between y/1/2 and 2yf§, and /x; is varied 
between the two natural choices, Mz and rtit, which is 
approximately a factor of two variation. A similar plot 
of the variation for ulul —> ZlH is shown in Fig.[3jThe 
only example with significant \ih variation is transverse 
gauge boson production. Fig. [4] shows the \ih variation 
for ulul — > W^W^ . The m dependence is still very 
small, but the Uh variation is nearly 10%. This is being 
investigated further. The high-scale matching is large for 
this process, almost 20%. We have checked that the fih 
dependence of the high-scale matching agrees with the 
EFT anomalous dimension. 

The scale dependence is below 1% in all cases except 
Wt production, so Uh and ui will be held fixed in subse- 
quent plots. The main contribution to the scale depen- 
dence is from the two-loop cusp anomalous dimension in 
the running, which is not cancelled by logarithms in the 
matching, since that has only been computed to one-loop. 
The scale-dependence would be even smaller if only the 
one-loop SCET anomalous dimension was used. The ui 
dependence increases with energy because of the log Q/ ui 
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FIG. 2: Variation of the ulul —> clCl cross-section as a 
function of Uh for [ii — Mz (solid red) and fn = mt (dashed 
blue), for y/s = 1 TeV (lower panel) and y/s = 5 TeV (upper 
panel). The cross-section has been normalized to that at y,h — 
y/s and hi — Mz- 




FIG. 3: Variation of the ulul — * ZlH cross-section. See 
Fig. [2] caption. 



FIG. 4: Variation of the ulul —> W^W T cross-section. See 
Fig. [2] caption. Note the tenfold change in vertical scale. 



term in the low-scale matching. The total radiative cor- 
rection increases even faster. 

The relative size of different contributions to the ra- 
diative corrections, such as the two and three-loop cusp 
and non-cusp anomalous dimensions, high-scale and low- 
scale matching, are shown for ulul — > clcl in Fig. [5l 
ulul -> W^Wt in Fig. u L u L -» W L W L in Fig.[7]and 
ulul — ► HZl in Fig. [8] We have determined the size 
of the electroweak corrections by turning them on and 
off in the effective theory, while keeping the full one-loop 
matching. 

In processes such as qq — > qq, the tree-level matching 
has terms of order on, i = 1,2,3 from B, W and gluon 
exchange, and the one-loop matching has terms of order 
ctiOij from box graphs. One has to decide whether to treat 
the ai 7 2 tree-level matching as an electroweak correction 
to the a3 term, and whether ai^ot^ terms in the one-loop 
matching should be treated as electroweak corrections to 
the tree-level QCD matching, or QCD corrections to the 
tree-level electroweak matching. To avoid this, we have 
used the full one-loop matching in all comparisons. The 
EFT graphs are one-loop diagrams, so it is clear whether 
to call them QCD or electroweak corrections. 

The electroweak radiative corrections are substantial, 
and increase rapidly with energy. For transverse W 
production, the electroweak corrections exceed 20% at 
Vs = 1 TeV. The corrections for transversely polarized 
gauge bosons are larger than those for longitudinally po- 
larized gauge bosons and the Higgs. At high energies, 
longitudinally polarized gauge bosons and Higgs bosons 
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FIG. 5: Different contributions to the UlUl — > clcl cross- 
sections as a percentage of the total rate. Lower panel: elec- 
troweak corrections (solid red), two-loop QCD cusp anoma- 
lous dimension (dashed blue), low-scale matching (dotted 
green), and the high-scale matching (dot-dashed cyan). Up- 
per panel: two-loop non-cusp QCD anomalous dimension 
(solid green), three- loop QCD cusp anomalous dimension 
(dotted blue), and two- loop electroweak cusp anomalous di- 
mension (dashed red). 



behave like members of the scalar doublet, and their 
electroweak radiative corrections depend on the SU{2) 
Casimir of the fundamental representation CV = 3/4, 
whereas transverse gauge bosons depend on the Casimir 
of the adjoint Ca = 2, which is larger. The high-scale 
matching for transverse W production is almost 20%, and 
both QCD and electroweak matching is important. The 
large \ih dependence is an indication that the two-loop 
high-scale matching might be important for this problem. 

The electroweak corrections are more important than 
the two-loop QCD cusp anomalous dimension. 3 The only 
corrections that are necessary for sub-1% accuracy are 
the one-loop QCD and electroweak corrections, the one- 
loop matching, and the two-loop QCD anomalous di- 



3 Note that a s /a 2 ~ 2.7 at 1 TeV and ~ 2.4 at 5 TeV. 



FIG. 6: Different contributions to the UlUl — * W^W T cross- 
sections as a percentage of the total rate. See Fig. [5] caption. 



mension. The two-loop electroweak and three-loop QCD 
corrections are completely negligible; even the two-loop 
QCD non-cusp anomalous dimension barely exceeds one 
percent in a few cases. The one-loop high scale matching 
can be as large as 10%. Higgs corrections proportional to 
the top-quark Yukawa coupling arc a few percent. These 
arise even for processes not involving top quarks via vir- 
tual top quark loops in the anomalous dimensions for the 
longitudinal gauge bosons and Higgs (see Ref. [HI]). 

Finally, in Figs. 151 RBI we show the complete radiative 
corrections to the scattering rate for some representative 
processes involving quarks, transverse and longitudinal 
gauge bosons, and the Higgs, normalized to the tree-level 
rate. At high energies, the radiative corrections suppress 
the scattering rate, but at lower energies, the radiative 
corrections lead to an enhancement. The rates arc shown 
as a function of vi for different values of t. The lower fig- 
ure shows the total correction to the rate, whereas the up- 
per panel shows only the electroweak contribution. The 
numerical values of the electroweak contribution are dif- 
ferent from the previous plots, because they have been 
normalized to the tree-level rate. The previous plots com- 
puted the electroweak contribution using the ratio of the 
total rate including electroweak corrections, to the total 
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FIG. 7: Different contributions to the ulul — > WlWl cross- 
sections as a percentage of the total rate. See Fig. [5] caption. 



FIG. 8: Different contributions to the ulul — > HZl cross- 
sections as a percentage of the total rate. See Fig. [5] caption. 



rate without elcctroweak corrections. 

Fig. Q3] shows the radiative corrections to gg — > 
WtWt, which is a background for the Higgs search at the 
LHC. This process first occurs at one-loop in the stan- 
dard model. We have not included the radiative correc- 
tions to the high scale matching C, which would involve 
computing two-loop diagrams, but all the other radiative 
corrections have been included using the effective theory. 



III. KINEMATICS AND NOTATION 

We will compute amplitudes for an arbitrary scatter- 
ing process with r external legs in the high energy limit. 
The large energy scale will be denoted by Q, and can 
be chosen to be the center of mass energy in the colli- 
sion. The effective theory computation is valid in the 
limit that the kinematic invariants pi -pj are of order Q 2 , 
and invariant masses are small, m, <C Q, Mw,z <C Q, 
where is the mass of particle i. The resulting ampli- 
tudes can be used, for example, to compute the rate for 
W+jets as long as the jet invariant masses are small com- 
pared with Q. The SCET computation neglects power 
corrections of the form M^ z /Q 2 an d m i/Q 2 in the one- 



loop and higher radiative corrections, and these are es- 
timated in Sec. [Vj The full dependence on mass ratios 
such as m 2 /M^, z , Mfj /M^ z , which are order unity in 
the power counting, is included in the EFT computation. 

It is convenient to introduce unit three- vectors rij, 
|n,-| = 1, which point along the direction of motion of 
the scattering particles. For incoming particles i, de- 
fine the null vectors rii = (1, Hi) and hi = (1,— ru), with 
rii-hi = 2. For outgoing particles j, rij and hj are defined 
by rij = — (l,iij) and hj = — (1, — n.j). 

In SCET, fast moving particles in the rii direction are 
referred to as rij-collinear particles, and are described by 
fields £ ni . A four- vector p can be decomposed as 



- (m + - (rii ■ p) nf + 



(1) 



If p is the momentum of an rii collinear particle, and 
p 2 = A 2 Q 2 , A < 1, then h ■ p ~ 0(Q), n-p~ C(A 2 Q), 
p± ~ O(XQ). If rii is chosen to point precisely along the 
direction ofpj, so that pj_ = 0, then hi ■ pt = Ei + \pi\ ~ 
2Ei and n t ■ p, = E t - |pj| ~ m 2 /(2E t ). 

In an r-leg scattering amplitude, all momenta pi are 
chosen to be incoming, so that for an outgoing particle 
j, Pj is the negative of the momentum of the particle. 
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FIG. 9: Plot of the ulUl —> clCl cross-section normalized 
to the tree-level rate, for t/s — —0.2 (dotted blue), —0.35 
(dashed red), —0.5 (solid black), —0.65 (dashed magenta) and 
—0.8 (dot-dashed cyan). The lower panel shows the rate with 
QCD and electroweak corrections, and the upper panel with 
only the electroweak corrections. 



FIG. 10: Plot of the ulUl — > W^W T cross-section normal- 
ized to the tree-level rate. See Fig. [5] caption. 



The labelling convention chosen for the Higgs doublet 



rij ■ pj ps 2Ej > is still positive for outgoing particles, 
because of our definition of rij for outgoing particles. 

A massive particle such as the top quark is described 
by a velocity four- vector v 1 *, with v ■ v = 1, where = 
7(l,/3), 7 = l/Vl — /3 ■ /3. For energetic top-quarks, 
it is sometimes convenient to use the four-vector (3^ = 
(l,/3), with (3 2 = I/7 2 — ► in the high-energy limit. 
This allows for a smooth transition in the high energy 
limit to a massless description, with [3 — > n. 

The Sudakov form factor will play an important role 
in this paper. The spacclikc Sudakov form factor F(Q 2 ) 
is defined as the particle scattering amplitude by an ex- 
ternal current, with momentum transfer Q 2 = —q 2 > 0. 
It is convenient to compute the form factor in the Breit 
frame (see Fig. [T5|) . where the particle is back-scattered, 
and the momentum transfer q has q° = 0. The Sudakov 
form factor is an r = 2 scattering amplitude, where the 
incoming and outgoing particle are identical. In the Breit 
frame, n\ = (l,n), ri2 = — (1, — n) so that fi\ = —n^ and 
n-2 = -n\. 



1 

V2 



tp 2 + itp 1 
v + H — iip 3 



(2) 



so that (p a cx iT a (</>) . The charged gauge and Goldstone 
bosons are 



1 



W ± = —= (W 1 T iW 2 ) 
v2 

+ t 1 - 2\ 

ip = -j= (tp T ) , 



(3) 



and the sign convention for the Z and photon fields is 



Z = cos 8 W W 3 - sin 9 W B . 
A = sin 6 W W 3 + cos 6 W B . 



(4) 



The SU(2) and U(l) fine structure constants are a.i and 
a.\ respectively, and the QED fine structure constant is 
a cm , with 



1 



1 1 

Q'2 Ot\ 



(5) 
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FIG. 11: Plot of the u L u L 
ized to the tree-level rate. See Fig. [5] caption 



W^W L cross-section normal- 



FIG. 12: Plot of the ulUl — > ZtZt cross-section normalized 
to the tree-level rate. See Fig. [5] caption. 



We will use the abbreviations 
aw 



sin tfw 



OL1 , 



sin 2 Ow cos 2 6w 



c w 


= cos 


Ow 


s w 


= sin 


Ow ! 


Lm 


= log 


M 2 

7^ 

Q 2 


Lq 


= log 


m 2 


Lq/m 


= log 


Q 2 

M 2 



(0) 



All couplings are in the MS scheme so that Eqs. (|5|6j) 
provide a definition of sin 2 Ow{n) in terms of ai^Qx). 

A generic gauge coupling is denoted by a, and the 
electromagnetic coupling constant is cv e m- 

Anomalous dimensions, form-factors and matching co- 
efficients will be expanded in powers of a/(4nr), 



m a 
7 = 7 4^ +7 



2 > (5) 



(7) 



IV. STRUCTURE OF THE LOGARITHMS 

The schematic structure of the radiative corrections is 

\ 



/ 1 
aL 2 



an 



oL 



a 2 L 4 a 2 L 3 a 2 L 2 a 2 L a 2 



3| 6 



a 3 L 



(8) 



where the first row is the tree- level result, the second row 
gives the one- loop corrections, etc., and all the terms 
are added to give the total amplitude. The standard 
terminology used in the literature is to call the first row 
leading order, the second row NLO, the third row NNLO, 
etc. and refer to the a n L 2n term as LL, the a n L 2n ~ 1 term 
as NLL, etc. Thus the a 2 L 2 term is referred to as the 
NNLL term at NNLO. 

In the leading-log regime, where aL ~ 1, the terms 
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FIG. 13: Plot of the ulul — > HZl cross-section normalized 
to the tree-level rate. See Fig. [5] caption. 



have the form 



/ 1 



m 



a 

1 a a 2 



(9) 



and terms of high order become important. Clearly fixed 
order perturbation theory breaks down completely when 
L get comparable to 1/a, and one needs to include all 
terms along and below the 45° diagonal in Eq. ©. 

The SCET formalism gives a result for the logarithm 



4 The loading double-log series exponentiates, and is convergent, 
i.e. e z has a series expansion with infinite radius of convergence. 



FIG. 14: Plot of the radiative corrections to the GG — > 
W^W^ cross-section. The lower panel shows the rate with 
QCD and electroweak corrections, and the upper panel with 
only the electroweak corrections. The radiative corrections 
are independent of t. 



v/VV\A/\/>- 



P2 
— Pi 



FIG. 15: The Breit frame, where the incoming particle is 
backscattered by an external current. 



of the amplitude in the form 

/ ctL 2 evL a 



log9Jt = 



or L a 
a 3 L 4 a 3 L 3 a 3 L 2 a 3 L a 3 
a 4 L 5 



(10) 



which is the known form of the perturbation series for the 
Sudakov form factor [(32|, [63|, [6J] . Eq. (fT0|) implies many 
non-trivial relations between the coefficients in Eq. |8|). 
In the leading-log regime, the exponentiated series be- 
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comes 



log m 



1 a 
1 a a 2 
1 a a 2 



a 



(ii) 



/ 



To get a reliable value for the amplitude in the leading- 
log regime, it is essential to sum the first two columns of 
Eq. (jlip . In this case, the higher order terms are small 
and under perturbative control, even though the total ra- 
diative correction can be very large (order I /a). We will 
use the exponentiated form of the amplitude. The first 
column in Eqs. (|9llip is the LL series, the second column 
the NLL series, etc. 

The SCET analysis automatically gives the amplitude 
in the exponentiated form Eq. (|10|) . By expanding the 
SCET result in a scries in a, one can recover the form 
Eq. ([8]). The exponentiated form must be used for the 
QCD corrections. It also needs to be used for electroweak 
corrections in W production at several TeV, where the 
electroweak corrections become order unity. At lower en- 
ergies (~ 1 — 2 TeV) it is not necessary to sum the entire 
set of terms along and below the 45° diagonal in Eq. j9]). 
It is still necessary to include the two-loop L 4 , L 3 , L 2 cor- 
rections, and in some cases, the three- loop L 6 , L 5 terms. 
The SCET formalism automatically sums all these terms 
(as well as higher order L enhanced terms). The two- 
loop a 2 (L 4 , L 3 , L 2 ) terms can be obtained using a one- 
loop renormalization-group improved computation. We 
have checked that the SCET formalism reproduced the 
known two-loop] L" n > 1 terms JnUB, El EI EI El, 

m Ei, Ei, m m m m, m, m, mS 

It is worth emphasizing that with the EFT theory 
method, no additional work needs to be done to sum the 
higher order terms. This is done automatically by the 
renormalization group. Instead one needs to work harder 
to drop these terms by expanding the SCET expressions 
in a, and then truncating the series. By contrast, in a 
fixed order computation, one first computes the one- and 
two-loop graphs, matches the resulting amplitude to the 
infrared evolution equation and then exponentiates the 
resulting amplitude, so summing the higher order terms 
requires additional work. 



V. ESTIMATE OF NEGLECTED POWER 
CORRECTIONS 

The EFT computation neglects 1/Q 2 power correc- 
tions in the one-loop and higher radiative corrections. 
There are also power corrections at tree-level, which arise 
from W, Z, m t masses in the propagators and particle 
phase space. These are easy to include, and do not limit 
the computation. Thus the neglected power corrections 
are small, because they only arise from loop graphs. 



To estimate the size of power corrections in the loop 
graphs, we compare the exact and EFT results for the 
one-loop Sudakov form factor with a massive gauge bo- 
son. The exact one-loop result is 



F(Q 2 /M : 
olCf 



1 



-in 
M 2 

Q 



7 2tt 2 4tt 2 M 2 , M 2 



M 2 , A 

Q 1 



3 Q 2 
M 2 



Q 2 



+ 2l °S-Ql l °S ( 1- 

M 2 , M 2 , 



M 2 

Q 2 " 



-2Li 9 



M 2 

Q 1 



M 2 



-Li 2 



M 2 

"q 2 " 



(12) 



In the Q 2 ^> M 2 limit, this becomes 



F(Q 2 /M 2 ) = F EFT (Q 2 /M 2 ) 
aC F M 2 



4tt Q 2 
« ( mA 

which is the EFT result. 



M 



M 2 



4_7T_ 2 

y 



(13) 



F EFT (Q 2 /M 2 ) 



= 1 



aCp 

47T 



2^ 2 M 2 



log 2 



Ml 
o 2 " 



(14) 



plus M 2 /Q 2 suppressed corrections. We define the power 
corrections 6F(Q 2 /M 2 ) as the difference F(Q 2 /M 2 ) - 
F EFT (Q 2 /M 2 ). At Q 2 = M 2 , 



SF(Q 2 = M 2 



uGf 

4-7T 



(15) 



In Fig. [TBI we have plotted the power corrections 
as a percentage of the trcc-lcvcl value F(Q 2 /M 2 ) = 
1, using an SU{2) gauge theory with coupling a = 
ow/(47rsin 2 6w), the value of the standard model SU(2) 
gauge coupling at Mz- The power corrections shown by 
the solid black line are less than 2% even at Q = M, and 
fall below 1% by the time Q/M > 2. The EFT radia- 
tive corrections, shown by the dotted red line, grow with 
energy. 

The other place where power corrections enter is 
through the use of the Goldstone boson equivalence theo- 
rem to compute the cross-section for longitudinally polar- 
ized gauge bosons, since the polarization vector is re- 
placed by the momentum k^. Here the power corrections 
exist even in the tree- level amplitude. In Fig. 1171 is plot- 
ted the ratio of the qq — > WlWl tree- level cross-section 
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FIG. 16: Plot of the one-loop power corrections to the Su- 
dakov form factor (solid black), the EFT one-loop correction 
(dotted red), and their sum, which is the total one- loop cor- 
rection (dashed blue) as a percentage of the total form factor. 
The gauge coupling constant has been chosen to be the stan- 
dard model weak SU (2) value. 



FIG. 18: Ratio of the qq — > WlWt cross-section to qq — > 
WlWl at tree- level, as a function of the center-of-mass energy 
for center-of-mass scattering angles of 45° (long-dashed blue) , 
60° (dotted red), 90° (solid black), 120° (short-dashed green) 
and 135° (dot-dashed brown). 



c 

0) 

u 
u 
a 



10 



■10 



-\ 1 \ 1 


i i i | i i i i | i i i i 


:\ \ 




\ 




- \ \ 








~ / 




T V 




- I: 
_ 1 
J! 

/',!,,,,! 


, , , 1 , , , , 1 , , , , " 



10 15 
E CM /(2M) 



20 



25 



FIG. 17: Percentage difference between the tree-level qq — > 
WlWl cross-section computed using the equivalence theo- 
rem and the tree-level cross-section computed directly, as a 
function of the center-of-mass energy for center-of-mass scat- 
tering angles of 45° (long-dashed blue), 60° (dotted red), 90° 
(solid black), 120° (short-dashed green) and 135° (dot-dashed 
brown) . 



AI 2 /El M relative to the qq -> W T W T and qq -> W L W L 
cross-sections. The WlWt cross-section has been omit- 
ted in our analysis. In Fig. [TH we plot the ratio of the 
WlWt to the WlWl cross-section. The WlWt rate falls 
with energy, and is sub-1% of the total cross-section at 
E CM /(2M) > 15. 

The parton-level results in this paper have an accuracy 
better than 1%, with the exception of Wt production. 
From the estimates in this section, we see that they can 
be used for processes other than those involving external 
longitudinal W and Z bosons at relatively low energies, 
2 — 3 times Mz, while retaining their sub-1% accuracy. 
For Wl, Zl processes, the use of the equivalence theorem 
introduces power corrections at tree-level, so the correc- 
tions are not sub-1% till much higher energies. One can 
achieve sub-1% accuracy at much lower energies by in- 
cluding the tree- level power corrections in the final result, 
and by including the tree- level rate for WlWt produc- 
tion. Then only a/ (An) suppressed power corrections are 
omitted, as for the Sudakov form factor, and the results 
can be used close to threshold with high accuracy. 



VI. THE EFFECTIVE THEORY AMPLITUDE 



ctet computed using the equivalence theorem, to the ex- 
act expression Uexact, as a function of the center-of-mass 
energy Ecm for different center-of-mass scattering angles 
#CM- The ratio CET/cexact asymptotes to one, but the 
corrections do not go below 1% until £'cm/(2M) > 15. 
The correction is ~ 8% for £'cm/(2M) = 5. Production 
threshold is B CM /(2M) = 1. 

The qq — > WlWt cross-section is suppressed by 



The next few sections contain a lengthy discussion of 
the general form of the scattering amplitude to all or- 
ders in perturbation theory. We have already seen that 
for practical purposes, a one-loop analysis augmented by 
the two-loop cusp anomalous dimension is sufficient. The 
general form of the amplitude to this order is known ex- 
plicitly, and the subtleties discussed in the subsequent 
sections are not relevant. We therefore provide a brief 
summary of the subsequent sections, and the form of the 
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amplitude. 

A generic r-particle scattering amplitude can have sev- 
eral different gauge invariant contributions. One writes 
the amplitudes and anomalous dimensions as matrices in 
a standard operator basis. For example, the qiq 2 — ► 9394 
amplitude can have the two possible color structures, 5 
written schematically as 



■03 ^4 



(16) 



where the color of the initial particles is combined to a 
color singlet, or color octet, respectively. The amplitude 
and anomalous dimensions are matrices in this operator 
basis. A more useful notation is to use operators in color 
space acting on the various particles [651 ] . rather than 
matrices in a particular basis. Tf is a color operator 
acting on particle i. The action of the color operators on 
the particle fields is 



A. The Renormalization Group Improved 
Amplitude 

At high energies, fixed order perturbation theory 
breaks down because of large logarithms of the form 
log(pi ■ pj/M 2 ), which are of order L ~ logQ 2 /M 2 in 
the log-counting, since all invariants pi ■ pj (i 7^ j) are 
of order Q 2 . From the point of view of the high-energy 
theory, these logarithms arise from infrared divergences 
as M — > 0. The same logarithms arise from ultraviolet 
divergences as Q — > 00 in the EFT, and can be summed 
using the EFT renormalization group equations. 

The renormalization group improved on-shell scat- 
tering amplitude (i.e. an 5-matrix element) has the 
schematic form [f| § 



(T A ) a ^ i0 8 ij = -(T A ^ i ) a 5 ij 



Tf^ a = fa (T A ) J l3 = {^T A )J K 



(17) 



where T A is the color matrix in the representation of 
particle i. The minus sign in the first of Eq. (|17p is be- 
cause ipi annihilates particle i. One can think of Tf as a 
field-theory operator Tjj acting by commutation, so that 



T A ,A, 



= - (T 



(18) 



Taking the hermitian conjugate then gives the second of 
Eq. (fT7|). Then, for example, 



Tf • Tf ^4 Mi 



-^3^4 ^T B T B ^ 



Tf • Tf ^T A ^ ^ 2 T A Vi 



-VS 3 T A V4 i> 2 T B T A T B ^ 



Cf - ^C A ) </> 3 T^4 VaT'tyi 



Thus Ti • T2 is equivalent to the matrix structure 



-C F 
- (C F - \C A ) 



(19) 



(20) 



in the basis Eq. (jl6|) . 

The number of color invariants grows rapidly with the 
number of external legs, and the color operator method 
provides a compact way of writing the group theoretic 
form of the amplitude. The soft anomalous dimension 
and matching has universal form when written using 
color operators. 



5 The term color refers to the SU (3) X SU (2) X U(l) gauge quantum 
numbers, not just to QCD. 



SPt = exp 



xPexp 



D ( a(/z ; ), log— ,* 



M 2 
7? 



df 1 



7(«(M).*) 



xc(a(^),{log^j) (21) 

where boldface quantities are matrices, and P denotes 
path ordering so that higher values of /j, are to the right. 
The Hh and fit scales are of order Q and M, respectively, 
and the \Xh,\ dependence cancels between the matching 
conditions and the renormalization evolution, so that StJl 
is fj,h,i independent. The invariants pi ■ pj are of order 
Q 2 . * denotes a possible dependence on kinematic vari- 
ables that is explained below. Further details on the form 
Eq. ((HD are given in Sees. IVlTTl H 

In Eq. (j2Tj) , c is the high scale matching coefficient 
at fih, and can be a function of all the invariants 
{log(p,; ■ Pj)/Hh}- It is independent of low-energy scales 
such as the gauge boson mass M or fermion masses 
m -C Q. The matching scale \Xh is chosen of order Q, 
so that {log(pi ■ Pj)/fJ-h} are all of order unity in the log- 
counting. 

The anomalous dimension 7(a(/x),*) is the SCET 
anomalous dimension used to evolve the operators from 
fih to in. 7 is given by the ultraviolet divergence struc- 
ture of the effective theory, and is independent of low- 
energy scales such as M and m. The * denotes linear 
dependence on log(pi - pj)//! 2 and a possible dependence 
on cross-ratios (also referred to as conformal ratios) of 
the form 



ijkl 



log 



(Pi ■ Pj){Pk ■ Pi) 
(Pi -Pk)(Pj - Pi) 



(22) 



which are order unity in the log-counting, can potentially 
arise in scattering amplitudes with four or more external 
particles beyond two-loop order ^SJBS Sim, HI, 
and are discussed further in Sec. |Xj 

In deep inelastic scattering and B — > X s j, the SCET 
anomalous dimension is linear in log Q 2 /fi 2 to all orders 
in perturbation theory [56l . |66| . and can be written as 
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7(m) = A («(//)) log/z 2 /<5 2 + B(a(/^)). A is known as 
the cusp anomalous dimension, and B will be referred 
to as the non-cusp anomalous dimension. The A term 
is order L in the log-counting, and is proportional to 
the unit matrix. The integral of this term is order L 2 , 
and gives the exponentiated Sudakov double-logarithms 
in high-energy scattering. The B term is order unity in 
log-counting, and its integral gives an order L contribu- 
tion to log 9Jt. B can have a non-trivial matrix structure 
in color space. 

The analysis in Sec. IVIIHXII is concerned with deter- 
mining the structure of the anomalous dimension and 
low-scale matching to all orders in perturbation theory, 
and determining what possible variables can occur in *. 
For the applications discussed in Sec. [Til it is sufficient to 
know the form of the anomalous dimension and matching 
to two-loop order. This is known by explicit computa- 
tion, and many of the complications discussed in later 
sections do not occur to two-loop order. We summarize 
here the results to two-loop order. The reader not inter- 
ested in the technical details of their derivation can skip 
directly from the end of this section to Sec. IXIIII 

The matching D and anomalous dimension 7 can both 
be written as the sum of collinear and soft contribu- 
tions, D = Del + D5 and 7 = 7cl + 7s, where the 
collinear contributions are unit matrices in color space. 
The collinear running and matching can be written as 
the sum of one-particle contributions, 

7c = ^JcA^ -Pi/v) > 

i 

D c = J2 D c,i(ni-Pi/^M/fx), (23) 

i 

where the sum is over the incoming and outgoing exter- 
nal particles in the hard scattering process. The collinear 
contributions are process-independent, so that, for exam- 
ple, 7c and Dc for a«i quark have the same values in 
ulUl — > dhdh and in u^dh — ► W + Z. The collinear 
anomalous dimensions have the form 

1C,i{ni-Pi/n) = r(a(/x))Ti -Tilog-^ + Bi(a(ji)) . 



/' 



(24) 



Recall that fn -pt is twice the energy of the particle. The 
first term is order L in the power counting, and the second 
term is order unity. T is the cusp anomalous dimension. 
At one- loop, 



r(aOO) 



(25) 



and the two-loop value is given by multiplying by the K- 
factor, Eq. ([55)1 . 6 T has no subscript i because it does 
not depend on the particle. 



The collinear-matching has the form 

TL ' ' T) ' 

Dc,i{fii-Pi/tj) = JiaiMlog-t-Z-Ti-Ti+Eiiain)), 



and at one-loop, 



J(a(/x)) 



2tt 



, M 2 
log — , 



(26) 



(27) 



and does not depend on the particle. The collinear 
matching at the low scale /i; has a dependence on the 
high scale Q through the log(fij -pi) term. J is related to 
the cusp anomalous dimension, and is required for con- 
sistency of the low-scale matching, see Sec. IXVI 

The matrix dependence is in the soft anomalous dimen- 
sion and matching. These have a universal form when 
written using the color operator notation, 



7s - -r(a(/i))^T I -T J log. 



-Hi ■ 71 j — I 



0+) 



m 

D s = -J(a(/i),L M )^Ti -T, log- 
to) 



(— n, ■ rij — i0 + ) 



(28) 



See Sec. IX Dl for the definition and references. 



where the sum is over all pairs of external particles in 
the scattering process. The two-loop values are given by 
multiplying by the A'-factor, as shown in Refs. [IH, [44( 
for the soft anomalous dimension in QCD. 

The soft-collinear split divides the amplitude into a 
collinear contribution which depends only on the ener- 
gies fii ■ Pi of the individual particles and is independent 
of their directions, and a soft contribtion which depends 
only on the vectors {rii} and hence on the directions of 
the external particles, but not on their energies. The 
full amplitude depends only on the Lorcntz invariant 
dot products pi ■ pj , and the dependence on the vectors 
{n-i} introduced by the effective theory cancels between 
the soft and collinear terms by reparamctrization invari- 
ance [63,111]. 

The anomalous dimension has a cusp contribution F, 
which is order L, and the rest, which will be referred to 
as the non-cusp anomalous dimension, given by Bi and 
75, collectively referred to as B. 

The low scale matching D = D^l + Do has Dl pro- 
portional to the unit matrix which contributes at order 
L, and Do with non-trivial matrix structure which con- 
tributes at order unity. Dl is given by the J term in 
the collinear matching, and Do by the Ei terms in the 
collinear matching, and the soft matching Dg. The low- 
scale matching does depend on the low-energy scales M 
and m. In the standard model, it depends on the elec- 
troweak symmetry breaking scale, and includes contribu- 
tions from the top-quark mass and 7— Z mixing. The new 
feature of the high-energy scattering problem first found 
in Refs. H,@,0] ^ s * ne presence of the single- log term Dl 
in the low-scale matching. This term can be shown to 
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Series 


r 


B 


-Dl Do c 


LL 


V 






NLL 


V 


V 


V V <V 


NNLL 


3^/* 


2^* 


2 V W 



TABLE I: The loop-order of various quantities needed to sum 
the LL, NLL, etc. series. F and B are the cusp and non-cusp 
anomalous dimensions, Dl and Do are the log and non-log 
parts of the low-scale matching, and c is the high-scale match- 
ing, y/ means the term has been included in our numerical 
results, y/ * means the term has been included, except for the 
scalar contribution. 



exist by explicit computation, and is required by consis- 
tency of the effective theory — the \ii dependence in the 
running from the cusp anomalous dimension T, which is 
of order L, is cancelled by the m dependence of D\. One 
can prove that the low scale matching is order L in the 

n theory, 
[ESQ. This is 
discussed further in Sec. IVIIII 

The leading-log (LL) series is given by the one- 
loop cusp anomalous dimension, the next-to-leading-log 
(NLL) series by the two-loop cusp anomalous dimension 
and the one-loop values of B and Di, the NNLL series 
by the three- loop cusp, two-loop B and D±, and one- loop 
Do and c, and the N n LL series by the n+1 loop cusp, 
the n-loop B and D-y, and the n— 1 loop Do and c. 

Eq. (f2"Tj) for the standard model sums the QCD and 
electroweak corrections, including cross terms such as 



log-counting to all orders in perturbation theory, so that 
terms such as log 2 Q 2 / (if do not appear 



a s a\ t 2, Oi s yf, or ai^yf which depend on mixed prod- 
ucts of the Yukawa, strong and electroweak coupling con- 
stants. Table |l] shows the order that the various quanti- 
ties are needed to sum the different series. The \J shows 
which terms have been included in our numerical results, 
and shows the terms which have been included ex- 
cept for the scalar contribution. The only quantities re- 
maining to sum the NNLL series are the log term in the 
low-scale matching Dl, and the scalar terms in the three- 
loop cusp and two-loop non-cusp anomalous dimension. 
The two- loop value of D L has been computed in Ref.|15| 
for fermions when mn = Mw, and is given in Eq. (|192[) . 



B. The UV/IR Correspondence 

The gauge theory under consideration is a sponta- 
neously broken SU (2) gauge theory with a Higgs in the 
fundamental representation, so that all the gauge bosons 
are massive. Our analysis is also valid in the more general 
case of an arbitrary broken gauge theory with massive 
gauge bosons. The advantage of studying this theory in- 
stead of a massless gauge theory such as QCD is that the 



infrared structure of the theory is under complete per- 
turbative control, with the divergences being regulated 
by the gauge boson mass M, if the gauge coupling at M 
is weak. 

The high energy behavior of a spontaneously broken 
theory is the same as that of the unbroken theory, since 
spontaneous symmetry breaking is soft, and can be ne- 
glected in the high energy limit. Thus our computa- 
tions also gives the high energy behavior of gauge theory 
amplitudes in an arbitrary unbroken gauge theory with 
massless gauge bosons. The ultraviolet behavior of an 
unbroken gauge theory is identical to the infrared behav- 
ior in perturbation theory, so we also obtain information 
on the infrared structure of gauge theory amplitudes in 
perturbation theory, which is a topic of current inter- 
est HE S3, m, m M EE Ell- 
in, the broken gauge theory infrared divergences are 
regulated by the gauge boson mass, so the renor- 
malized amplitudes are functions of a, Lq and I_m, 
A(a(jj,), Lq, I_m), and the scattering amplitudes are IR 
finite. In a massless gauge theory, the scattering am- 
plitudes are infrared divergent, and can be regulated in 
perturbation theory using dimensional rcgularization, so 
they have the form A(a(n), Lq, e) with 1/e playing the 
role of I_m- The form Eq. (|2"T]) holds for the massless 
gauge theory, with the low-scale matching replaced by a 
low-energy parton matrix element computed in pertur- 
bation theory. This is the same as Eq. (|21[) with I_m — > e 
in the D terms: 



9H = exp 



xPexp 



T>(a(fj,i),e,* 
^ d^ 



Pi ■ Pj 



7 («(//),*) 



xc ( a{/j, h ), <J log 2 



(29) 



where 7, and c are unchanged from their values in the 
broken gauge theory. The 1/e terms in Eq. (|29[) are in- 
frared divergences, and arise because one is computing 
on-shell scattering amplitudes with a fixed number of ex- 
ternal particles. They cancel between real and virtual 
graphs if one computes experimentally measurable quan- 
tities such as jet production rates, which are infrared 
safe. 

We will illustrate the relation between the massive am- 
plitude Eq. (|2"Tj) and the massless amplitude Eq. (|2"9")) in 
Sec. IXVI by deriving the Magnea-Sterman equations for 
the Sudakov form factor for both theories, and com- 
paring the results. 



VII. FACTORIZATION REGULATOR 

In SCET, the collinear and soft integrals need an in- 
frared regulator, which cancels in the total amplitude. A 
regulator is necessary for all factorization calculations — 
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otherwise one can show that the logQ 2 //z 2 term in the 
anomalous dimension, which leads to the Sudakov double 
logarithms, would be absent. 

A convenient regulator is the A-regulator introduced 
in Ref. [4l[. The A-regulator for particle i is given by 
replacing the propagator denominators by 



1 



1 



(Pi + k) 2 - to 2 (pi + k) 2 - mf 



Aj has mass dimension 2. This regulator can be im- 
plemented at the level of the Lagrangian, since it corre- 
sponds to a shift in the particle mass. The on-shell con- 



production. We will use the boosted HQET (bHQET) 
formalism of Refs [H| , in which heavy particles are 
described by HQET fields W^hvj, where vi is the ve- 
locity four-vector Vj = 1 and nj is a null vector n 2 = 
along the direction of particle I. 8 The momentum is 
Pi = mivi, and 7/ = £j/toj = nj ■ uj/2 is the boost. 
W ni contains ultracollinear fields that couple to h Vl anal- 
(30) ogous to the collinear Wilson line W ni for collinear fields 
£rii • W ni does not require a regulator after zero-bin sub- 
traction [4ll ] . 

For the h v propagator and vertex, Eq. (|31I32[) become 



dition remains p~ 



The form Eq. (f5U|) determines 



the regulator to be used for soft and collinear graphs in 
the effective theory. 

In SCET, the rij-colhnear propagator denominators 
have the replacement Eq. (|3H)) . Collinear gauge invari- 
ance requires a collinear Wilson line W ni , so that rii 
collinear fields only appear in the combination [W^.£ ni ] . 
As shown in Ref. [U, [6^|, this is true even in the 
presence of the A-rcgulator, since the collinear Wilson 
lines do not require a regulator after zero-bin subtrac- 
tion [U|6^, [zO, |7l|j[23- 7 raj-collinear graphs only depend 
on the regulator Aj, and are independent of Aj for j =/= i. 

If particle i interactions with a soft-gluon (fc M — > 0), 
the interaction with a A-regulator becomes 



2tf 



2tf 



(31) 



(pi + k) 2 - A, 2pi -k-Ai' 
Using Eq. ([T]), and keeping only the leading term gives 
(fit ■ pi) n? 



(rii ■ pi) (rii ■ k) - Ai (th ■ k) - S t 



where 



rii - Pi 



(32) 



(33) 



It is convenient to use Si as the regulator for the soft 
gluon interactions, even though it is determined from Aj 
by Eq. (f3"3"]) . Si has mass dimension 1. 

The soft interaction Eq. (|32[) is invariant under the 
rescaling rii Si^ii, Si — ► SiUi, so that soft graphs can 
only depend on the combination nf/Si. This can also be 
seen directly from 



nf/Si 



(«! ■ k) — Si {rii/ Si ■ k) - 1 



(34) 



so that only the ratio /Si appears in soft graphs. 

We also need the A-regulated form for graphs contain- 
ing heavy particles to study processes such as top quark 



2rf 



(pi + k)' 2 



A, 



vi-k- A' 



where 



Aj_ 

2mj ' 



(35) 



(36) 



and 



so that the HQET propagator is regulated by A'j, 
amplitude only depends on the ratio vi/A'j. 

The soft interactions of the heavy quark are given in 
terms of a soft Wilson line Y ni , for which Eq. (|31I32|) 



become 



(nj -pi)n^ 



{ni ■ pi) (ni ■ k) - Ai (raj ■ k) - Si ' 



where 



A, 



A, 



2A'j 

hi ■ pi mi(hi ■ vi) hi ■ pi 



(37) 



(38) 



and the amplitude only depends on the ratio rij/Si. A'j 
and Si for heavy quarks both have mass dimension 1. 

In previous papers, we have used the analytic regula- 
tor [2J,[2l, defined by letting 



(39) 



p 2 — m 2 



(p? - TO 2 ) 



1+Si 



where Si arc dimcnsionless. The same regulator is also 
used in Feynman graph computations using the method 
of regions (see, e.g. Ref. [HI)- The analytic regulator is 
scaleless, so that SCET soft-graphs vanish on-shell. How- 
ever, the analytic regulator breaks the color Ward iden- 
tities necessary to write the collinear fields in the form 
[Wn£n] @, if there is more than one gauge invari- 
ant amplitude. It can used for the Sudakov form factor 
(where there is only one color invariant), while maintain- 
ing the [W^£„] form of the collinear fields. 



7 The necessity of subtractions has been shown previously using 
other regulators |73|| . 



We will use lower case subscripts i for light particles and upper 
case subscripts / for heavy particles. 
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VIII. FACTORIZATION STRUCTURE OF THE 
AMPLITUDE 

The scattering amplitude Wl computed in the effective 
theory is computed from soft and collincar diagrams. 971 
will refer to the amplitude, or to the anomalous dimen- 
sion, depending on the context. We will use rm to denote 
the masses and any other quantum numbers of particle 
i (such as whether it is a fermion or scalar), and to 
denote its gauge quantum numbers. The argument {T 1 / 4 } 
denotes dependence on the gauge quantum numbers of all 
the external particles, whereas denotes dependence 
on the gauge quantum numbers only of particle i, etc. 
Note that all amplitudes depend on the total particle 
content of the gauge theory, for example through the (3- 
function and vacuum polarization diagrams. What {Tf} 
and T/ 1 denote is dependence on the quantum numbers of 
the particular external particlc(s) for a given scattering 
amplitude. 

In the general r-particle scattering amplitude, all the 
momenta and gauge indices will be treated as flowing into 
the vertex. Thus an outgoing particle with momentum 
p in representation D\ is treated as an incoming particle 
with momentum — p in representation 91. The analysis 
below holds for an operator at arbitrary momentum, so 
we do not assume that the total incoming momentum 
J2iPi vanishes, e.g. for the Sudakov form factor (r = 2), 
momentum transfer q is inserted at the vertex. 

A generic r-particle operator has the form 

r 

f({h})\[{Y ni ) bi a > [Wlt ni ] ai (40) 
i=l 

where a* are gauge labels, and / is a Clebsch-Gordan 
coefficient to make a gauge singlet operator. Here £ n; is 
an n.;-collinear external line which couples to ni-collinear 
gauge fields, but not to soft gauge fields. W ni is the rii- 
collincar Wilson line in the representation of particle i 
necessary for rii-collinear gauge invariance, and Y ni is a 
soft Wilson line in the n, direction in the representation 
of particle i and contains soft gauge fields. The structure 
of the interaction is shown in Fig. \T§[ and shows how the 
gauge indices are contracted. Soft gluons couple only to 
the Y n s, the triple lines in Fig. \T§[ and they can interact 
with any of the triple lines, as shown in Fig. 1201 rij- 
collincar gluons can interact with W ni and £ ni , so they 
interact with the double line and dashed line for a given 
direction rij. They do not communicate between two dif- 
ferent directions. 

The rii collinear interactions involve only the [W^ ^ ni ] 
part of the operator. Thus they only depend on the in- 
teractions of particle i, and are independent of all the 
other particles in the process. Note that the Wilson line 
W ni is in the gauge representation of particle i, and does 
not need a A-regulator. Thus collincar gauge invariance 
implies that each collinear amplitude has the form (to all 
orders) 

cxpli (aM.LMJ/'.A.y,™^)!, (41) 




FIG. 19: Structure of the r-particle operator in the effective 
theory. The dashed lines are collinear fields £„, the double 
lines are collinear Wilson lines WX, and the triple lines are 
ultrasoft Wilson lines Y n . The gauge indices on the r ultrasoft 
Wilson lines are combined into a color singlet at the vertex. 




FIG. 20: Soft gluons can interact with all the Y n s (triple 
lines). 



where is an abbreviation for all the quantum num- 
bers of particle i, such as its mass, and whether it is a 
fermion, scalar, or gauge boson. The amplitude is a unit 
matrix in gauge space since [Wn<£*i<] gets corrections 
which do not change the index . The index cannot mix 
with the other indices, since the other indices do 
not transform under n.;-collinear gauge transformations. 
The exponential form is convenient for the subsequent 
analysis. The collinear amplitude depends on the A reg- 




/ S 

/ s 
/ \ 
/ \ 
/ \ 
/ \ 



FIG. 21: n-collinear interactions couple W n (double line) and 
£n (dashed line). 
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FIG. 22: The soft function can depend on m t through vacuum 
polarization graphs. 

ulator only for particle i, which enters into the particle-? 
collinear propagator. 

The soft amplitude is given by graphs in which the 
gauge bosons couple to the soft Wilson lines Y ni via 
eikonal interactions. The soft amplitude can be written 
(to all orders) as 

cxpS (a(At), L M , {Tf} , {/< A}) . (42) 

The amplitude is a matrix in gauge space, since the soft 
interactions can mix the gauge labels in Eq. (|4"0|) . so exp S 
is a matrix exponential. The amplitude is written as 
exp S rather than S for future convenience. S depends on 
the gauge quantum numbers and directions n., of all the 
particles. The argument in Sec. I VIII shows that it only 
depends on n,; and Si in the combination n^/Si. Note 
that the soft graphs do not depend on the masses of the 
external particles, or whether they are fermions, scalars, 
or gauge bosons. Thus the amplitude does not depend 
on {rrii}. The soft interactions do not generate e a 0\o-, 
and so can only depend on the invariants (n, • Tij)/(Si8j). 

The soft amplitude can depend on particle masses 
through vacuum polarization graphs such as those in 
Fig. 1221 They arise from the mass-modes discussed in 
Ref . [57], [H| . This dependence is omitted from the argu- 
ment of Eq. (|4"2"|) for simplicity, and does not change the 
following discussion. The dependence on particle masses 
{m,} is an invariant dependence, i.e. it does not depend 
on the particular external legs in the scattering ampli- 
tude. 

It is well-known that the Sudakov form factor depends 
on Q 2 = 2pi-p2, and that the r-particle scattering ampli- 
tude depends on log pi -pj. This can be seen explicitly by 
a one-loop computation. If the soft and collinear sectors 
were truly independent, then the scattering amplitude 
would be completely independent of pi ■ pj . This is be- 
cause each collinear sector is Lorentz invariant and can 
only depend on one momentum at a time, whereas the 
soft sector does not depend on {pi}, so there is no way 
to generate pt ■ pj . The flaw in the argument is that the 
collinear and soft sectors need a regulator to be sepa- 
rately defined, and communicate via the regulator. The 
Pi ■ pj dependence is generated through this regulator de- 



pendence. This is not specific to SCET, but is a univer- 
sal feature present in any method which factors the total 
amplitude into soft and collinear contributions. All such 
factorization methods must necessarily introduce a fac- 
torization scale that acts as a regulator and enters both 
the soft and collinear functions. The pi ■ pj dependence 
enters through this regulator dependence, as we will see 
explicitly below for the A regulator. 
The total amplitude is 

r 

971 = exp S JJexpftl] , (43) 

i=l 

and the matrix ordering between Ii and S is irrelevant 
since the collinear Ii terms are proportional to 1. The 
amplitude is the product of the collinear and soft contri- 
butions since the different sectors are independent. The 
log of the scattering amplitude is 

log 931 = S^( M ),L M ,{^},{^^|) 

r 

(44) 

where matrix ordering between collinear and soft is again 
irrelevant because the collinear terms are proportional to 

a. 

The total amplitude is independent of the infrared reg- 
ulator, since the SCET scattering amplitude in a massive 
gauge theory is not infrared divergent. This is because 
the effective theory has the same infrared behavior as 
the full theory, and the full theory has no infrared diver- 
gences since they are all regulated by the gauge boson 
mass. This implies that the regulator dependence must 
cancel between the soft and collinear terms. 9 The only 
reason for the A-regulator is that divergences are intro- 
duced in the effective theory by the split into collinear 
and soft sectors. 

The cancellation of the regulator in Eq. ((44)) puts 
strong constraints on the form of the effective theory am- 
plitude. The S dependence must cancel between the soft 
function and the r collinear functions. The collinear am- 
plitudes are a sum of one-particle terms, each depending 
on a single Si , whereas the soft amplitude depends on the 
products SiSj of all possible pairs of Ss. 

The soft function depends on r(r — l)/2 variables (rij • 
rij)/(SiSj). [Remember that nf = 0.] It is convenient to 
rewrite these in terms of r(r — 3)/2 cross-ratios which are 



9 One can write down the integrals for the soft and collinear am- 
plitudes, and first combine the integrands before evaluating the 
integral. This allows one to compute the total amplitude with- 
out the necessity of introducing a regulator such as the A regu- 
lator miii. 
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independent of {Si} and r variables which depend on Si. 
Define 



C(ij\r-2,r- l,r) = log 



(rij ■ rij)(n r —2 ■ n r )(n r _i • n r ) 
{rii ■ n r ){rij ■ n r ){n r - 2 • n r -\) 

(45) 



for 1 < i < j < r and 

C(ir\r — 2, r — 1, r) 







(46) 



for i < r. In Eqs. (|45l46p . the particles r — 2, r — 1 and 
r play a special role, and have been included explicitly 
in the indices of C. Most of the time, we will keep these 
reference particles fixed, and denote the cross-ratios by 
C(ij). C(r — 2, r — 1) = 0, so there are r(r — 3)/2 non- 
trivial independent cross-ratios. The cross-ratios Eq. (|4"5"|) 
are homogeneous in {rii} and so have the same values if 
rii — ► 7ii/&. 

The cross-ratios are invariant under rescaling rii ~ * 
XiUi, and so one can make the replacement ni — ► (rii ' 
Pi)rii, i.e. rii — > p^, and write them in terms of ratios of 
dot products of momenta, 

C(ii|r-2,r-l,r) - log ft ' gfcj ' ' . 

CP* • Pr ) (Pj ■ Pr ) (Pr-2 • Pr-l) 

(47) 

The cross-ratios first appear for r = 4, and so can enter 
into the 2 — > 2 scattering amplitude. They are absent 
for the Sudkaov form factor (r = 2). For two-particle 
scattering 1 + 2^3 + 4, the only cross-ratios are 

7>i 2 43 =C(12|234) = log 

, \yx fi)\fA ft) , 
= lo § 71 TT7I — = lo S ■ 

^1342 =C(13|234) = log 

= lo § 7 \1 \ = lo § ~o 



(ni 


• " 2 )(n 3 


■ n 4 ) 


(m 


■ rii){m 


■ n 3 ) 


(pi 


P2)(P3 ■ 


Pa) 


{pi 


Pi)(P2 ■ 


P3) 


(ni 


■ m)(n 2 


■ ni) 


(ni 


■ ni){n 2 


■ n 3 ) 


(Pi 




Pa) 


(Pi 


Pi)(P2 ■ 


Ps) 



(48) 

in terms of the Mandelstam variables. 
Define the remaining r variables by 

y(i\r -2,r-l,r) = ^^gj r { s 

2 (n r -2 • n r )(n r -i ■ n r ) 

+ log(ni • n r ), i < r 

1 n r~2 ' n r-l 

y(r\r-2,r-l,r) = --log- r. 

2 (n r -2 ■ n r )[n r -i ■ n r ) 

(49) 

In terms of these variables 

Iog(»i • n 3 -) = C(ij) + y(i) + y(j), i ? j (50) 



if we let C(ij) = C(ji) if i > j. Thus any cross-ratio of 
the {rii} which is invariant under rescalings of the rii can 
be written in terms of C: 



log 



(rij ■ nj)(n k ■ n{) 
(rii ■ n k )(nj ■ m) 



C(ij)+C(kl)-C(ik)-C(jl) 

(51) 



since y(i) cancel on substituting Eq. (|50|) . 

Then the r(r — l)/2 variables \og(fi 2 rii ■ rij)/(SiSj) can 
be written as a linear combination of C(ij) and y(i) — 
\og(Si/fi), so the soft function S can be considered to be 
a function of C(ij) and X(i) = y(i) - log^/p). 10 The 
collincar function Ii is a function of A^//i 2 = Si (rii -pi)/ ^ 2 
(see Eq. (J33j) ) , so it can be regarded as a function of 
Z(i) = log(5i/fi) + \og(hi - pi/fx). The total amplitude 
can then be written as 

r 

logOJt = S({C(ij), #(*)}) +5^(201 (52) 

i=i 

where only the momentum and regulator dependence has 
been shown explicitly. Since VUt is independent of Si, 
differentiating it w.r.t. Si gives 



as 



dX(i) dZ(i) 



1 



(53) 



Since only depends on Z(i), differentiating w.r.t. 5j 
gives 







d 2 S 



dX(i)dX(j) 



(54) 



and all mixed partial derivatives of S vanish. This implies 
that S must have the form 



S = J^Si(X(i)), 



(55) 



a sum of matrix functions each depending on a single 
X(i) variable. Substituting Eq. (|55|) in Eq. (|53")l gives 



dX(i) dZ(i) 



(56) 



from which it follows that both sides are a constant since 
X(i) and Z(i) can be varied independently. The constant 
must be proportional to the unit matrix since the r.h.s. 
is. It can only depend on variables other than X(i) and 
Z(i) which are common to both S and Ii and have been 
suppressed in Eq. ([56]) . Thus one gets 



dX(i) dZ(i) 



1 = J (a((i),L M ,TA) 1 (57) 



It is useful to use Si/fi so that the variables are dimensionlcss. 
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so that J only depends on the color of particle i and the 
gauge boson mass, and is independent of particle masses, 
and whether they arc scaiars, fcrmions, or gauge bosons. 
Thus J is a universal function independent of i, since the 
color dependence is explicit through the argument T A . 
Integrating Eq. ([57]) gives 



h = 



and 



J (a(fx),L M ,T A ) lX(i) 
+G l (a(p),L M: {T t A } ,{C(i 3 )}) , 

J (a(n), \- M ,T z A ) lZ(i) + E (a(n), L M ,T A , rm/n) , 

(58) 



log Tt = ^j(a(M),L M ,lf)l 



y(i) + 



n r ■ Pi 



+ ^S(a(^),L M ,i; A ,m l //i) 1 

i=l 

+G («(/*), L M , {Tf},{C(ij)}) , 



(59) 



where G = J^G*. The {Si} dependence has cancelled, 
and the amplitude depends linearly on the combination 
X(i)+Z(i) =y(i)+log(n l -p i ). The final result Eq. (59]) 
is completely symmetric under the exchange of particle 
labels. This is not manifest in Eq. ((59|) because the split 
of kinematic variables into C(ij) and y(i) is not symmet- 
ric. The final expressions we need for our computation 
given in Sec. IXIVl are manifestly symmetric. 

The form Eq. |59|) is valid to all orders in perturbation 
theory. J is related to the cusp anomalous dimension. 
The factorized form Eq. (|59j) is valid for the renormal- 
ized amplitude logStTt and for the anomalous dimension 
7 in the effective theory, both of which have the addi- 
tive structure Eq. ([52)) . The EFT anomalous dimension 
is given by the UV divergences in the EFT, and does not 
depend on IR scales such as M and to.; , so these variables 
can be dropped for -y. The anomalous dimension can de- 
pend on fii -pi and on the directions n$, since these appear 
explicitly as labels on the EFT operator. The amplitude 
Eq. (|59|) is linear in log (i to all orders in perturbation 
theory, a result first shown in SCET for deep-inelastic 
scattering in Ref. [56|. For the multiparticle case, the 
log \x dependence is proportional to the unit matrix, and 
commutes with the other contributions. 



A. The Light-Light Sudakov Form Factor 

The above discussion works for four or more external 
legs, r > 4, where one can use the variables C(ij\r — 2, r — 
l,r) and y(i\r - 2,r - l,r) defined by Eqs. (|45|49|) . For 
r = 3, there are no cross-ratios, but one can still use 
Eq. (|49|) and Eq. (f50|) . dropping C(ij) in the final result 
Eq. p|. 

For r = 2, the soft function depends only on the sin- 
gle variable (ni ■ H2)/(#i#2)j and there is only a single 



color invariant, so that there is no matrix structure. 
The derivation leading to Eq. (|63|) for the two-particle 
case has already been given earlier in several different 
ways [!, H, @, HI]. We repeat the derivation because 
we will use the results for r = 2 (the Sudakov form 
factor) extensively in the remainder of this article. In 
Eq. |52|) , the soft function can be chosen to be a function 
of X = log(/! 2 ni ■ n2)/ (28182)- Since 93T is independent 
of 8i t 2, differentiating w.r.t. 81,2 gives 



dS d/W 



dX dZ(l) dZ(2) 
Differentiating again w.r.t. 81^2 gives 
d 2 S 



dX 2 



0. 



(60) 



(61) 



so that 

S 
h 



J(a(fj,), L M , T t A )X + G(«( M ), L M , T t A ) , 
J(a(/x), L M ,T A )Z(i) + E(a(n), L M , T A 



m 



(62) 



Note that because the operator is a gauge singlet, Cf = 
Ti • T 1 = T 2 • T 2 = -Ti • T 2 and T A + T A = 0, so 
that there is only one independent color operator. We 
will drop the subscript i on T A when we write the color 
dependence of the Sudakov amplitude. 

Eq. (|60[) implies J\ = J2, so there is no subscript i on J. 
This also follows from the discussion below Eq. (f56|) . The 
collincar function Ii has the same value whether particle 
i scatters via a two-particle operator, or via an r-particle 
operator, so J and E are the same as in Eq. (|58| . 

Adding the pieces gives 



log 9JT = j(a(n),L M ,T A )log 



+ ^ J B(a(^),L M ,T A ,TO l / M ) 

+G(a(n),L M ,T A ) . (63) 
In deriving Eq. (|63p we have used J\ = J 2 = J and 



J 



log 



J log 



ni ■ n 2 
2 

-Q 2 
m 2 



log 



^1 - Pi 
/i 



log 



(64) 



where — Q 2 



{P1+P2) 2 = («i-Pi)(n 2 -_p 2 )(ni-n 2 )/2, 



since all momenta are incoming. 



B. Change of Basis Transformation 



The basis for cross-ratios Eq. (j45|) picked out particles 
r — 2, r — 1 and r. One can use a different reference set 
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of particles for the cross-ratios, r — > sq, r — 1 — > si and 
r — 2 — > s 2 in Eq. (|45|) . The cross-ratios in the new basis 
are given in terms of the old basis by 

C(ij\s sis 2 ) = C(ij\r - 2, r - 1, r) 

+C(siso\r — 2,r — 1, r) — C(isi|r — 2,r — 1, r) 
-C(js |r -2,r- l,r) , 



^(i 7^ s |sosiS2) = y(i\r - 2, r - 1, r) 

+ 2<?(s2Si|r — 2,r — l,r) - -C(s 2 sok — 2,r - l,r) 



can be written as 



1 



C(sis |r - 2, r - 1, r) + C{is \r - 2, r - 1, r) , 



3 ; (so|soSis 2 ) = -^C(s 2 si|r- 2,r- l,r) 

+ 2 C ( s 2So|r - 2,r— l,r) + ^C(sis |r - 2,r - l,r) . 

(65) 

The transformation is linear in y, and C does not involve 
y, so linearity of Eq. (|59|) in 3^ is maintained, and chang- 
ing the basis does not provide any additional constraints. 
If G is linear in {C(ij)} in one basis, then it remains linear 
in {C(ij)} in any other basis. Linearity in C is connected 
with Casimir scaling, and will be important in Sec. |Xj 



IX. THE SCATTERING AMPLITUDE 



7 ( M ) = a^r(«( M ),T/ 



i=l 



y(i) + log 



rij ■ Pi 



-l£f2(a(/« A , S i) 
-E (aW,{T/},{C(ij)}) 



(67) 



r, ft, E are the analogs oi J, E and G for the anomalous 
dimension. Since 7 does not depend on low scales such 
as M or rrii, these variables are dropped in Eq. (|6"T)) . T, 
fii and E are obtained from the l/e ultraviolet divergent 
terms in Eq. (|59[) . and obey the factorization structure 
of the EFT amplitude. T, which is the cusp anomalous 
dimension, was shown above to only depend on the color 
of particle i, but not on whether it is a fcrmion, gauge 
boson, or scalar. Thus we have written it as r(a(/Lt), T A ) 
where the color dependence is explicitly through the sec- 
ond argument. The function T is then universal, i.e. there 
is no need to introduce different functions Tj for each par- 
ticle. While f2 does not depend on the mass of particle 
i, it can depend on whether the particle is a fcrmion or 
scalar or gauge boson, so we have included an argument 
Si (for statistics) in f2; again there is no need for different 
functions fli. 

For the Sudakov problem, the scattering of particle i 
by an external current, r = 2, there is no color matrix 
structure, and 



One can now write down the general expression for a 
hard scattering amplitude using the EFT method. 

(1) One matches from the full theory on to SCET at 
a scale fx^ of order Q. This gives a column vector of 
coefficients c of the allowed SCET operators. We will 
normalize c so that it is dimcnsionlcss. c is a column 
vector, and has the form 



c(a(nh),{pi -Pj},Vh) 



(66) 



and can depend on the invariants pi ■ pj which are order 
Q 2 , and logarithmically on fih, but not on low scales para- 
metrically smaller than Q, such as M . Since c is dimen- 
sionless, it can depend on ratios such as (pi -Pj)/(pk 'Pi); 
which are invariant under an overall rescaling pi — > Xpi 
of the momenta, but are not invariant under individual 
rescalings pi — > XiPi. The invariance under pi — > Xpi 
is broken by the Hh dependence, i.e. by the anomalous 
dimension of c. c is computed from the hard part of 
the diagram, i.e. by using pure dimensional regulariza- 
tion with all infrared scales set to zero. This is the usual 
procedure for computing a matching condition, and the 
EFT factorization methods do not apply to this part of 
the computation, fx^ is chosen to be of order Q to mini- 
mize the logarithms in c. 

(2) Run the coefficients from ph down to a low scale pi 
of order M using the SCET anomalous dimension. The 
SCET anomalous dimension has the form Eq. (|59p . and 



0=2)/ \ 



r(a(M),T^) 
+a(a(p),{T l A }) 



1 & 
l0 V 



i=l 

(68) 



using the results of Sec. IVIII Al and the form Eq. (|63[) . 
E for the Sudakov problem, which we require later, has 
been called a, and the log(— 1)T term has been absorbed 
into a. The branch cuts are determined by analytically 
continuing the form factor from Euclidean to Minkowski 
space, or by an explicit graphical computation. While 
there is no color matrix structure, the anomalous dimen- 
sions till depend on the color of the particles (through 
CV), which is why the T, dependence is retained in 

Eq. (HHD. 

The evolution of the amplitude from fi^ — > pi is given 

by 



H 



Pexp 



d/i 



70) 



(69) 



where P denotes /i-ordcring with larger values of pL on 
the right. 

(3) Compute the low-scale matching at (i = yU; ~ M. 
The matching is also constrained by the EFT form 
Eq. (59]), 



D 



log d 
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l5>(a(/i),L M ,lf) 



y{i) + log 



n t ■ Pi 



+ H^£(a(^),L M ,!7f,™ i A«) 

i=l 

+G{a( f i),L M ,{T l A } 1 {C(ij)}) 



(70) 



D is linear in logQ (the log(ni ■ pi) terms) to all orders 
in perturbation theory. The existence of this logarithm 
in the low-scale matching,and linearity to all orders, was 
first pointed out in Ref. H, Q . D and d are square ma- 
trices. Previous studies of factorization have all been 
for masslcss gauge theories. Eq. ((70)) is the factoriza- 
tion form for the low-scale matching in a massive gauge 
theory, which has not been studied previously. 
We break D into the single- log term and the rest, 



D 

D L 



= Do 



= Xy(a(M),LM,7f)log 



i=l 



TH ■ P i 



D = l^J {a(n),L M ,Tf)y(i) 

r 

+lJ2E(a(^),L M ,T l A ,m l /fi) 

i=i 

+G{a(n),L M ,{T l A },{C(ij)}) , 

so the low-scale matching can be written as 

d = cxpD cxp(Z?il) 

where the two terms commute. 

For the Sudakov form-factor (r = 2) for particle i, 



(71) 



(72) 



Q2 

A),i + Dn log — 2 



A 

D lfi = j(a{fi) 7 L M: T t A ) , 

2 

D 0ti = ^ E (a(fi), Lm, T a , mj/jj) 

i=l 

+g(a([i) : L M ,{T A }) , (73) 

where G for the Sudakov problem is denoted by g, the 
log( — 1) J term has been absorbed into g, and there is no 
matrix structure. 

The final result for the amplitude is (showing the scales 
at which the various pieces are computed) 



971 = d(/x ; ) H(pi «- fx h ) c(n h ) 



(74) 



and is independent of \ih and \xi . This is the amplitude 
given in Eq. (|2"Tj) . 

It is worth emphasizing that the total amplitude fJJt in 
the gauge theory does not have the factorization proper- 
tics studied here. The factorization properties only hold 
for the effective theory amplitude, and hence for d and 



H in Eq. (|74|) but not for the hard matching correction 
c. Thus the factorization properties are obscured in an 
explicit diagrammatic computation. The effective theory 
reorganizes the computation to make the factorization 
structure manifest, analogous to the way in which heavy 
quark effective theory (HQET) makes heavy quark spin- 
flavor symmetry manifest [78[. 



A. Consistency Conditions 

Hh independence of Eq. (|T4[) gives 

Mt^ = 7(m)c(m) 
d/x 



and /i; independence gives 
eld 



so 



e _D >— e Do + tfi- 



(75) 



(76) 



-7(A0 ■ (77) 



d/it d/i 

Comparing the log(nj • pi) terms on both sides gives 

dJ , , 

>6j ^ T (78) 

which relates the \og(fii - pi) term in the low-scale match- 
ing to the cusp anomalous dimension. Thus the presence 
of the log(ni ■ p{) term in the low-scale matching follows 
from consistency of the effective theory Q . 



X. CASIMIR SCALING 

The results of the previous section simplify dramati- 
cally at one-loop. The collinear graphs are proportional 
to the color factor Tj • Tj and the soft graphs are a sum 
over pairs of two-particle terms proportional to Tj • . 
The collinear and soft amplitude can be written as 

Ii = Tj • Tj [/(a(//), Lm, Aj//i 2 , TOj/^)] 



(v) 



SiSj 



(79) 



where (ij) is a sum over pairs of particles with i ^ j. This 
is a special case of Eq. (|44| , where the color dependence 
of the collinear and soft amplitudes has been explicitly 
written, so that Ii and S do not depend on color, and 
S has no matrix structure. Furthermore, the S function 
multiplying Tj ■ Tj can only depend on the single variable 
(rii ■ n,j)/(5i5j), rather than all possible dot products, 
and so the EFT amplitude cannot depend on cross-ratios, 
which involve at least four particles. 

The EFT anomalous dimension and matching 
Eqs. (|103I104I105I106|) simplify considerably because of 
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the restricted form Eq. (|79p . It is simplest to rcderive the 
anomalous dimensions and matching conditions which 
follow from the ^-cancellation conditions starting from 
Eq. ([75]) , rather than imposing the restriction of Eq. ([79]) 
on the results of the previous section. The answer is 



log 971 = J(q(ai), L m ) J2 Ti ' Tl lo S 



-J(a(ji), l M )J2 T * ' T i l °S — 5 



E(a(n), L M , mi/fj,) + -g(a(fi), L M ] 



(80) 

and the derivation is given in Appendix The functions 
on the r.h.s. no longer depend on the color quantum num- 
bers, and the matrix function G has been replaced by a 
single function g. 



A. Anomalous Dimension 

Assuming Casimir scaling, Eq. (|75|) . the anomalous di- 
mension Eq. (|67[) simplifies to 

a>i - Pi 



7 (a0 = W^Ti.Ti log 



/' 



i0 + 



(U) 



(81) 



and the Sudakov anomalous dimension Eq. (|68[) becomes 
7 f =2) (M) = ^ • Ti 



+fl (a(n),s 2 ) + a(a(n)) 



(82) 



using Q 2 = (m • Pi)(n 2 • P2)(-ni • n 2 )/2. 

It is convenient to break the anomalous dimension into 
a collinear and soft anomalous dimension, the precise def- 
inition of which is given in Sec. IXII for the general case. 
The collinear anomalous dimension is 



7c 



r(t*(A0)io { 



rii - pi 



(83) 



the soft anomalous dimension is 



ls = -r(a( M ))^T,-T,log 

(ij) 



(—Tii ■ rij — i0 + ) 



(84) 



and the total anomalous dimension is 7 = 7c + 7s- The 
soft anomalous dimension for an r-leg amplitude has a 
universal form, and is completely determined in terms of 
the cusp anomalous dimension. At one-loop, 



Km) 



(85) 



Eq. |84|) is also valid at two-loops. 



B. Low-Scale Matching D 

Assuming Eq. (|75|) . the low scale matching Eq. ([70)) 
simplifies to 



D = logd 



= J2 J(a(M),LM)Ti-Ti log 



i=i 



rij ■ Pi 



-J(a( M ),L M )E T >- T >g- 



rii ' rij — 1 



0+ 



r 

E T *' T * E(a(p), L M ,mi/(i) 
i=l 

l^Ti-Ti g(a(ii),L M ). 



(86) 



D is linear in logQ (the log(n,; - pi) terms) and in log/i 
to all orders in perturbation theory. 

We break D into the single-log term and the rest, so 
Eq. (j7"Tj) becomes 



rii - pi 



D = D + 1£>l, 



D L = J2 J(a(ji),L M )Ti -Ti log 



D = -J(a(^),L M )E T * - T J l °S- 



E T *' T * E(a(fi), 

i=l 

■i^T.-T, g(a(fi),l M ), 



so the low-scale matching can be written as 

d = cxpD cxp(DlI) 

where the two terms commute. 
For the Sudakov problem 

O 2 

A = Am + -Dialog 

M 



(87) 



(88) 



A,i 
D t 



J(a(/i),L M )T i -Ti 



Lm, mi /fi) 



-9 {a(n), L M ) 



(89) 
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The matching can also be written as the sum of a 
collinear contribution 



J(a(n), L M ) lo{ 



+E(a(p), L M ,mi//i) + -g{a{ji), L M ) 



(90) 



and a soft contribution 
D s = 



JMmUm^TVTjIos 



(— rti • rij — i0 + ) 



(91) 



with D = T>c + Dg. TTie so/t matching for an r-leg 
amplitude has a universal form, and is completely deter- 
mined in terms of the single-log term in the matching, 
which is related to the cusp anomalous dimension. At 
one- loop, 



J(a(/i), U 



ftp") I 
2tt 



2-7T jti^ 



(92) 



The above expressions give a remarkable simplifica- 
tion in computing scattering amplitudes. The collinear 
contributions are a sum over one-particle terms which 
do not depend on the scattering process, and the soft 
contributions are universal. The soft functions have the 
minimal form necessary to restore reparametrization in- 
variance [6?], [H| (i.e. boost invariance) of the EFT am- 
plitudes. For back-to-back particles n.j ■ nj — —2, and 
log(-TH • nj)/2 = 0. 



C. Form of the EFT Functions 

Comparing the formulae in this section with 
Eqs. (|67l68l71l89p . we see that 



E (a(fi), l M ,Tf,mi/fi) 
9i{a{n), \-M,Tf) 



T, 
Ti 
T ; 
T, 
T, 



T 4 r(a(/i)) 

T 4 a(a(ji)) 

T t J(a(n),l M ) 

Tj E(a(n), Lm, rrii/fj,) 

T ifl (a(/i),L M ) (93) 



so that the color dependence is factored out into a pref- 
actor. In addition, the non-trivial matrix structures S 
and G in the soft contribution are determined in terms 
of scalar functions a and g which appear in the Sudakov 
problem. 



D. Higher Order 

The one-loop functions Eq. and the anomalous 
dimension Eq. ([52")) and amplitude Eq. (|8"6")) have two 




FIG. 23: Three loop correction to the cusp anomalous dimen- 



very special properties. The first is proportionality to 
Ti ■ Tj, which is referred to as Casimir scaling. Obvi- 
ously, this is a trivial result at one-loop, since the one- 
loop graphs all explicitly have this form, but it is non- 
trivial at higher orders. The cusp anomalous dimension 
T(a(fj,), T/ 1 ) has been computed to three-loop order [79j . 
and obeys Casimir scaling to this order. The cusp anoma- 
lous dimension is 

T(a(^Tf) = £(£fr(«)(T^), 



r« = 

r (2) = 
r (3) = 

where the ^-factors are 
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The scalar contributions to K\ is known [3, UB] , but not 
to K 2 . 

Casimir scaling can be seen from the proportionality 
of r^ 2,3 - 1 to r^ 1 -*, so that all three terms in T are propor- 
tional to C F . It appears that the C F term in K 2 violates 
Casimir scaling; however, it originates from the two-loop 
vacuum polarization correction Fig. 1231 and C F t F n F is 
from Tr (T a T fc T c T b ), and so the same for all particles, 
and does not depend on the particular particle whose 
cusp anomalous dimension is being computed. 

The non-cusp part of the collinear anomalous dimen- 
sion B = 20, + a does not obey Casimir scaling. The one 
and two-loop terms for fermions are [lj, [l5| 
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and the C|. term in B^ violates Casimir scaling already 
at two loops. In color operator notation, it is propor- 
tional to (Tj • Tj) 2 . The three-loop value B^ is known 
for a gauge theory with only fermionic matter [7!| . 

The cusp anomalous dimension is the anomalous di- 
mension of a Wilson loop with a point where the tangent 
vector is discontinuous (the cusp) [80| . Casimir scaling 
of Wilson line expectation values has been studied exten- 
sively in the literature, and there has been recent interest 
in this question 0, 0, |H, M, S, 0, M, ED, ■ Casimir 
scaling of Wilson loops has also been studied on the lat- 
tice, to test whether the string tension is proportional to 
the color Casimir of the quark. Recent results seem to 
indicate there is deviation from Casimir scaling [8lT ] for 
the string tension. Casimir scaling of Wilson loops does 
not hold when non-perturbative effects are included. For 
example, large Wilson loops in the fundamental represen- 
tation are confining and have an area law (in the absence 
of light fcrmions), whereas those in the adjoint represen- 
tation are screened, and have a perimeter law. In the 
large N limit, one can relate the two, since screening is 
1/N 2 suppressed [HHl]. 

The collinear anomalous dimension respects Casimir 
scaling at least till three loop order, by explicit computa- 
tion [73] . Recently, Aybat et al. [43|, 0] have shown that 
the soft anomalous dimension matrix at two-loop order 
also respects Casimir scaling, i.e. 



(2) 

Ts 



(97) 



where K\ is the same factor as for the collinear anoma- 
lous dimension, and is a number. Thus the entire two- 
loop anomalous dimension matrix can be obtained by 
multiplying the one-loop matrix by (1 + K\a/ {A-k)) 1 as 
can the Dl term in the two-loop matching. 

For practical LHC computations, we only need the one- 
loop electroweak corrections, the one-loop QCD match- 
ing, and the two-loop QCD running. The one-loop terms 
have the nice Tj • T, and sum-on-pairs form. The results 
of Aybat et al. [HI, 0] tell us that this form can also be 
used for the two-loop anomalous dimension. The two- 
loop non-cusp collinear anomalous dimension does not 
obey Casimir scaling, but it is a one-particle term and 
can also be trivially included in the results. 

An interesting question currently being investigated 
is whether Casimir scaling continues to hold at higher 
orders. Becher and Neubert [5l[ have shown that the 
soft anomalous dimension should continue to have the 
Tj • Tj\og(pi ■ pj) form at three- loops, and have sug- 
gested that this form persists to all orders. If this form 
holds at higher orders, then the form derived in Sec. IXII 
in terms of collinear and universal soft functions will con- 



tinue to hold at higher order. For the electroweak prob- 
lem, we need not only the anomalous dimension, but also 
the low-scale matching D, and it would be interesting to 
investigate the form of D at higher orders. The Becher- 
Neubert [5l[ argument that there are no color structures 
at two and three-loops other than Tj ■ Tj implies that 
the low-scale matching should continue to have this form 
to three-loops. 



E. 



Sum on Pairs 



The results in Sees. IX AIX Bl have two properties - 
they are proportional to Tj • Tj, and they are linear in 
log rii ■ rij ; , or equivalcntly in log pi ■ pj . Linearity in log pi ■ 
Pj implies that the amplitude can only depend linearly on 
the cross-ratios C(ij); higher powers would induce (log Pi ■ 
Pj)(logpk - pi) terms. The entire EFT one-loop result, 
including the low-scale matching, has linearity in both 
Tj • Tj and in logpj • Pj. The color structure follows 
trivially from the color factors of a one-loop graph, and 
linearity in log pi-p_j can be seen by explicit computation. 
As shown in Ref. pj , an amplitude of this form can also 
be written as a sum-on-pairs over two-particle Sudakov 
form factors weighted with Tj • Tj , 



log9Jt 



(ij) 



TjF(- Pi - P j). 



(98) 



Thus all the information in the scattering amplitude is 
contained in the Sudakov form factor. 

Becher and Neubert [5l| give a proposed form for the 
all orders anomalous dimension, which is both linear in 
Tj • Tj and in log pi ■ Pj, so that the sum-on-pairs for- 
mula Eq. (f9"8"|) would hold to all orders. Ref. Q used the 
analytic regulator for SCET. On-shell soft graphs vanish 
with the analytic regulator, so the entire amplitude arises 
from collinear graphs. The analytic regulator has the big 
disadvantage that the n-collinear sector fields cannot be 
combined into the W^£, n form of a single Wilson-line; in- 
stead one has separate n-collinear Wilson lines for each 
particle not in the n-direction, and the different collinear 
sectors have interlocking color structures 4l| . This com- 
plication does not hold for amplitudes with a single color 
structure, such as the Sudakov form factor, or particu- 
lar scattering processes such as 3 + 3 — > 6 + 8 in QCD 
(quix production). There is also effectively only a sin- 
gle color structure in the large- AT limit, because color 
reordering of the operators is suppressed by 1/N. The 
mixing is only suppressed for a suitable basis of opera- 
tors, those written using products of color matrices in 
cyclic order. 11 Thus in the qq — > gg example considered 



11 In the double-line notation [84j . each color connected part of the 
operator should consist of quark lines with alternating directions 
for the arrows, which are color-paired as one moves around the 
vertex in cyclic order, so that one has a planar vertex [83, 85]. 
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in Ref. [H], one should use the basis qq5 ab , qT a T b q and 
qT b T a q rather than qqS ab , qT c q f abc and qT c qd abc for 
the operator basis. One can then see that the anomalous 
dimension matrix in Ref. [86| becomes diagonal in the 
large- N limit. In Ref. 0], it was shown that for ampli- 
tudes with a single color structure, linearity in logp; • pj 
holds to all orders. Thus a check of the Becher-Neubert 
proposal requires a four-loop computation away from the 
TV — > oo limit. 

Linearity in log pi ■ pj implies that the amplitude and 
anomalous dimension have the form 

Am n i ' Pi 



log Tl = tJ2 </(a(M),L M ,lf )log 



where 



+lJ2E(a( f j,),l M ,T % A ,m i / f j l ) 

i 

-^Bfa) W ,L M )log - ni ^- iQ+ 
+G(a(/i),L M ) (99) 



2J2~B UJ) + J(a(n),L M ,lf)l = (100) 



and do not depend on the cross-ratios. Imposing in ad- 
dition linearity in T; • Tj for the logn^ • rij terms and 
in Tj • Tj for the rest reduces the amplitude back to 
Eq. ([90]) . Only using linearity in T ?; • Tj gives the general 
form Eq. ({59")) with a restriction on the color structure 
of the terms, but no constraint on the dependence on 
cross-ratios. 



XI. COLLINEAR AND SOFT FUNCTIONS 

The factorization form Eq. (|59[) of the EFT scattering 
amplitude allows us to efficiently compute the radiative 
corrections to high energy scattering processes, such as 
those needed for the LHC. The factorization structure 
studied in Sec. lVIlJ provides a more compact way of writ- 
ing our previous EFT results, by dividing the EFT scat- 
tering amplitude into a sum of one-particle collincar con- 
tributions, and a universal soft function. The procedure 
we follow is that used by Aybat, Dixon and Sterman [44[ 
to separate the amplitude into jet and soft functions. We 
will use the name collinear function, rather than jet func- 
tion, since the term jet function has already been used in 
the SCET literature for a different quantity. This decom- 
position will be used both for the anomalous dimension 
and for the low-scale matching D. 

The SCET computation automatically breaks the am- 
plitude into a product of collinear and soft contribu- 
tions. The splitting, however, is regulator dependent, 
and each piece depends, for example, on <5j if one uses 
the A- regulator. This is a generic feature of any factor- 
ization procedure which cannot be avoided, because the 



logQ 2 term in the SCET anomalous dimension is gener- 
ated precisely by this regulator dependence. It is nicer 
to reorganize the amplitude into regulator-independent 
terms which have the same factorization properties as the 
SCET collinear and soft contributions. The basic idea is 
to write the effective theory amplitude as the product of 
the square-roots of the Sudakov form factor in SCET for 
each external particle, and the rest. The square-root of 
the Sudakov form factor in the effective theory (i.e. the 
form factor without the high-scale matching) is called 
the collinear contribution, and the rest is the soft con- 
tribution. Since the total amplitude and the SCET Su- 
dakov form factor are both regulator- independent, this 
provides a regulator independent split of the amplitude 
into collincar and soft parts, which, as we will see, re- 
spects the factorization structure of Eq. (|59"1) . Our def- 
inition of the collinear function is identical to that of 
Aybat et al. [44[ for the jet function. Aybat et al. de- 
fine the soft function in terms of a ratio of web functions 
regulated by taking the directions n.i away from the null 
direction, and then taking the light-like limit. We have 
used a different regulator and defined the soft function 
in terms of SCET amplitudes, but the final result is the 
same. 

The Sudakov form factor for particle i is the amplitude 
for p\ — > p2, where both particles are of type i, and 
fi\ ■ p\ = t%2 ■ P2 = n% • Pi in the Breit frame. 12 The 
collinear anomalous dimension is defined as one-half the 
anomalous dimension for the amplitude, Eq. (|67p . and 
can be written using Eq. as 

7C?; = r(a( M ),Tf)log(^0 +QH U ),T/ l , Si ) 



/' 
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Half the log of the Sudakov low-scale matching gives the 
log of the collinear matching 



D 



= jAy+jA.ilog— ^j— 
= J( a ( M ),L M ,^)log^i 



+E (a(/z), \-M,T A ,mi/n) 



(102) 



using Eq. (ggj) . 

Now consider an arbitrary hard scattering process with 
r external lines, where the particles have energies 2Ei ~ 
n.j ■ pi . We will define the collincar anomalous dimension 
and collincar low-scale matching by summing Eq. (jlOip 



12 i.e. the Sudakov form factor is evaluated at the same energy that 
particle i has in the r-particle scattering process. An alternate 
choice is to evaluate all the functions at some common momen- 
tum transfer, such as 
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and Eq. ()102[) evaluated at • pi , and summed over all 
external lines. They are proportional to the unit matrix 
in color space. The total collincar anomalous dimension 
is the sum of one-particle terms 



7c 



r(a(/x),^)log 



+ -a{a(ii),Tf) 



as is the collincar low-scale matching 



(103) 



J(a(/i),L M ,I^)log 



-£7 (a(ju), LMjlV^jmi/zi) 



+ 25(«(m)> L M,T t ) 



(104) 



The soft anomalous dimension 75 is defined as the dif- 
ference of the total anomalous dimension and collincar 
anomalous dimension. Using Eq. (|67|) , we find 

r 

7 sM = l£r(a(M),35 A )l r (i) 

i=l 

-il£>(a( M ),^) 

i=l 

+£(«(//), {2^} ,{C(ij)}) . (105) 

The log of the soft matching D5 is the difference of the 
logs of the total matching D and the collinear matching 
D c , 

r 

D s = l^j(a(/i),L Ml 7^)y(i) 
»=i 



-l^ 5 (a(/i),L M ,{i; A },{C(ii)}) 

i=l 

+G («(/*), L M ,{i; A },{C(y)}) (106) 



and the total anomalous dimension and log of the low- 
scale matching are 



7 = 7c + Is , 
D = Dc + Ds 



(107) 



The decomposition of the anomalous dimension and 
matching into soft and collinear functions discussed 
above is that used by Aybat et al. [44| • The soft / collinear 
split is related to the decomposition of the SCET ampli- 
tude into soft and collinear graphs, but is not identical. 
The basic difference is that in an r-leg process, the soft 
contribution of the Sudakov form factors for each exter- 
nal particle has been added to the collincar function and 
subtracted from the soft functions (the a and g terms in 




FIG. 24: Structure of the r-particle operator in the effective 
theory for external light and heavy quarks. The dashed lines 
are collinear fields £ n , the solid lines are HQET fields h v , 
the double lines are collinear Wilson fines WX, and the triple 
lines are ultrasoft Wilson lines Y n . The gauge indices on the 
r ultrasoft Wilson lines are combined into a color singlet at 
the vertex. 



Eq. (|103I104I1051 HDSjl ). This cancels the regulator de- 
pendence of the SCET collinear and soft functions, and 
provides a regulator-independent definition for the two 
functions. 

The collinear functions are process-independent, and 
can be computed directly from the Sudakov form factor. 
The cusp part T of 7c and J of Dc are of order L in log- 
counting, whereas the remaining terms of 7c, Dc, as well 
as 7s, are order unity. One can move the log(nj -pi) 
terms from the collinear factors into the soft factors by 
rewriting n, • rij —> (pt ■ pj)/ fj 2 . It is more convenient to 
leave the order L terms in the collinear functions, where 
they are explicitly proportional to the unit matrix. In 
addition, the rij • rij form for the soft function remains 
valid for heavy quarks, so that the same soft function can 
be used for SCET and HQET fields. 

Assuming the Casimir scaling form Eq. (|79[) allows 7s 
and Ds to be written in a much simpler form. Using the 
general form given in Sec.|X]f or the anomalous dimension 
and amplitude, the collinear and soft anomalous dimen- 
sions and matching reduce to Eqs. (|80I84I90I85[) given in 
Sec. [Xj This is the form we need for the NLL standard 
model computations, which use one- loop matching and 
two-loop running, for which Casimir scaling is known to 
hold. 



XII. HEAVY QUARKS 

We can also compute the scattering amplitude for r + s 
external legs, where r legs are light, and treated using 
SCET fields W^Ui, != l,...,rands > 1 legs are 
heavy and treated using boosted HQET fields (bHQET) 
Wl T h vl , I = 1, . . . ,s [EE HI. The SCET operators arc 
given by Fig. 1241 where, for a heavy quark, the external 
dashed line is replaced by a HQET field. 

The rii collincar graphs have the same form as before, 
Eq. (141 [) . The on-shcll HQET graphs are also a sum over 
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the external heavy particles, 

expf// a(n), Lm,T/ 



A 'i ) 



1. 



(108) 



The HQET graph Hj does not depend on the heavy parti- 
cle mass, since it does not enter the HQET computation. 
The soft amplitude is given by graphs in which the gauge 
bosons couple to the soft Wilson lines Y ni and Y ni via 
cikonal interactions. The soft amplitude can be written 
(to all orders) as 

cxp S (a( M ), L M , {T a } , {/mf A} , {^/Sj}) . (109) 

The two amplitudes only depend on the combinations Vj, 
n^/Si and n^/Si as shown in Sec. IVIIl 
The log of the scattering amplitude is 

logSttt 

, Ax (f^m-rij n 2 rii-ni ^ 2 n I -n J 



8,5 j 



SiSj 



(a( f i),l M ,T t A ,A l /fi 2 ,m l / t i) 1 

i=l 
s 



(110) 



i=i 



where matrix ordering is again irrelevant because the 
collinear and HQET terms are proportional to 1, and we 
have written the n and v dependence in terms of Lorentz 
invariant dot products. 

By a now familiar argument, the cancellation of the 
regulator in Eq. (|110[) puts constraints on the form of the 
effective theory amplitude. Since the form Eq. (| 1 1 0[) is 
identical to Eq. (|44j) when there are no massive particles, 
the results of Sec. IVIlil hold. The only difference is that 
for heavy quarks, the collinear functions are replaced by 
HQET functions 

Hi = j(a( (J ,),L M ,Tf)tZ(I)+E h (a( l i),L M ,T A ) 

(111) 

analogous to Eq. (|58|) , where Z(I) = logA' 7 //j = 
log(£j//x) + log(n/ • vj)/2. Eh does not depend on toj, 
since mi is not a parameter in bHQET. Eh is not the 
same function as for a massless quark, but J is. Eq. (|59"|) 
becomes 



where hi-pi — > hi-vi/2 = 7/ because the relation between 
A ^ and 5i is now Eq. . The soft function G depends 
on color factors or cross-ratios involving all the particles, 
and any three particles (heavy or light) can be chosen as 
the reference indices in defining y and C. 

Equation Q112p is the general form for the SCET ampli- 
tude and anomalous dimension, and can be used for pro- 
cesses involving any combination of heavy and light par- 
ticles. The change on replacing a light quark by a heavy 
quark is n ■ p — > n ■ v and E — > Eh in the collinear func- 
tion. Thus the results of Sees . 1X1 and IXII can also be used 
for heavy quarks with this substitution. 



A. Matching from SCET to bHQET 

In processes where Q>m> M, the external particle 
is treated as an SCET field for Q > fj, > m, and as a 
bHQET field for m > \i > M, and we need to compute 
the matching correction at the scale fi = fi m ~ m be- 
tween the two theories. This result is need, for example, 
to compute QCD corrections to tt production. 

The difference between graphs in the theory above and 
below [i m is in the treatment of n-collinear interactions 
involving the heavy particle. In the theory above fj, m , 
they are given by W^ n , and in the theory below fi m , 
by W\h v . In the amplitude Eq. (fTTDj) . the differ ence is 
that Ii is replaced by the corresponding Hi for each ex- 
ternal heavy particle. Thus for each particle making the 
transition from SCET to bHQET, the matching factor is 



Vi 



Ii (a(fi),e,T z A , Ai/^i 2 -> 0,mi/fi) 



(113) 



logSJJl = £j(a(/«),L M ,2^)l 



y(i) + log 



n-i - pi 



nj ■ V] 



^j(«( Ai ),L M ,T / A )l^(/)+log v 

1 = 1 

r 

^2E(a0i),L M ,T i A ,m i /n)l 

i=l 
s 



1=1 



G (a(jx),L M ,{T A },{C(ij)}) 



(112) 



The computation in the theory above fj, m is done on-shell 
with p 2 = m 2 , and since m; 3> M one can set M = to 
compute the matching. The A-regulator is not needed, 
and dimensional regularization is sufficient to regulate 
the infrared divergences. The difference Ii—Hi is infrared 
finite. 

The matching Vi only depends on the collinear func- 
tions for light and heavy particles. The soft function 
remains unchanged at the SCET — > bHQET transition. 
The one-loop matching conditions are given in Table HT1 
and were obtained previously in Ref. [Q|. Two-loop re- 
sults for i-quarks can be found in Ref. [87| . 



B. The Heavy-Light Sudakov Form Factor 

The heavy-light Sudakov form factor is the scattering 
of an incoming HQET field (particle 1) to an outgoing 
collinear field (particle 2) by an external current, similar 
to the b — > S7 decay process studied in Ref. [lj. The 
result is given using Eq. (|112[) for the case r = s = 1. 
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47T 
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L m + £ + 2] 


B x 


a rp 
4tt A 
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"il 2 - 
2 Lm 


L m + £ + 2] 



D. Form Factor Relations 



TABLE II: The SCET -> bHQET matching correction for 
fermions, scalars, and transverse gauge bosons. L m = 
logm 2 /^ 2 , where m is the particle mass. 



The anomalous dimension is 



7 0*) = T(a(n),T A ) log 



(n x -pi){n 2 • v 2 ) 



+fi (a(/i), T A , si) + ft (a(ji), T A , s 2 ) 
+a(a( f i),T A ) , (114) 

where S1.2 also encode whether the particle is heavy, and 
we have dropped the subscript on T since 
Note that all heavy particles are the same, i.e. it does not 
matter if they are fermions, scalar, or gauge bosons. The 
low-scale matching is 



D = A, + £>i log 



(ni -pi)(n 2 • v 2 ) 



£>i = J {a(fi),l Ml T A ) , 

.Do = E h (a{^),L M7 T A ) +E(a(n) 7 L M ,T A 7 m 2 /^) 

+. 9 (a(/i),L M ,{T A }) . (115) 

C. The Heavy- Heavy Sudakov Form Factor 

The Sudakov form factor for two heavy external lines 

is 

7 (/i) = T(a( t i),T A )log(2v 1 -v 2 ) 

+Q (a(ti),T A , 8l ) + Q (a(fi),T A , s 2 ) 

+a (a( M ),T A ) , (116) 

and the low-scale matrix element is 

D = D + r>ilog(2«i • v 2 ) , 
Di = J (a(ix),L M ,T A ) , 
D = E h (a(fi),L M ,T A )+E h (a( l i),L M ,T A ) 

+g{a(n),L M: {T A }) . (117) 

The heavy-light anomalous dimension and matching 
are the average of the heavy-heavy and light-light values. 



From the results in Sees. IXII B[ IXH CI we see that 
Fh 1 ^h 2 (Q 2 )Fii^i 2 (Q 2 ) 



(118) 



where ?i,2 and ft, 1.2 are light and heavy particles, and 
h.2, hi >2 have the same color. This relation is satisfied by 
the EFT form factors. It is also satisfied by the full form 
factors in the limit Q 3> mi,2, M, and 



Em-i — >ni2 (Q ^)Erri3 — tin^iQ ) 



Ejni — '7714 (,Q 2 *)F n 



AQ 2 



' = 1 + 0^1,(119) 



because the high scale matching condition is independent 
of particle masses, and cancels out, as does the matching 
between SCET and HQET for each external heavy line. 
Eq. (|119p holds regardless of whether Q, m and M are 
widely separated, Q ^> m ^> M, or m and M are com- 
parable. So it also holds for the massless gauge theory 
M — > with the infrared divergences cancelling in the 
ratio. 



XIII. THE EQUIVALENCE THEOREM 

The radiative corrections at high energy for longitudi- 
nally polarized massive gauge bosons can be most easily 
com put ed using the Goldstone boson equivalence theo- 
rem [H, [M IM IM US HI US , which relates the gauge 
boson amplitude to that of the corresponding (eaten) 
Goldstone boson, 



In 



(A2...A?F\S\A» 
{-i£) n (i£) m {v ai . . 
' M 
~E, 



...A b L -I) 



+0 ^ 



(120) 



up to corrections which vanish in the high energy limit. 
The l.h.s. is the ^-matrix for an initial state to scatter to 
a final state, and the r.h.s. is the corresponding ampli- 
tude (made more precise later) with the gauge bosons re- 
placed by the corresponding (unphysical) Goldstone bo- 
son. £ is a radiative correction factor which is one at 
tree-level, and Ei are the energies of the gauge bosons. 
The existence of £ was first pointed out by Bagger and 
Schmidt [38|. The equivalence theorem only holds for 
the ^-matrix for physical states. All longitudinal gauge 
bosons A a h must be simultaneously replaced by the corre- 
sponding ip a for the theorem to be valid. Thus I and F 
do not contain any longitudinal gauge bosons. 

We briefly review the derivation of the equivalence the- 
orem including the radiative correction £ , and then give 
the result for £ in the MS-schcme. 

For a longitudinal gauge boson in the n direction with 
energy E, the momentum is p M = (E,pn) with E 2 = 
p 2 + M 2 hys , and the longitudinal polarization is chosen 
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to be 



FIG. 25: Two-point function. 



(p,£n)/M physj e L -p = 0, e| = -1, where 



•Wphys is the physical mass of the gauge boson. In the 
high energy limit e£ p»/M phys + 0(M/E). 

The general form of the two-point function in the Gold- 
stone boson/longitudinal gauge boson sector is 



T L P L 



ik v T^ A 



r T P T 



//// 



-ik ll r A *> 



(121) 



where the projection operators are 

> ) 



I' 1 ' 
T 



k 2 

9»v - 



(122) 



The momentum k flows into the blob through the line cor- 
responding to the right index of the matrix (see Fig. 
At tree-level, 



r T 


= M 2 - 


-k 2 , 






r L 


= M 2 - 


-h 2 

C ' 








= Y tpA 


= m(i- 


?)■ 






= k 2 - 


^' )2 M 2 , 


(123) 



where £ and £' are the parameters in the gauge-fixing 
term, and £ = £' in i?£ gauge at tree-level, and we have 
used the notation of Rcf. [40|. The tree-level Goldstonc 
boson mass is M£' / ■ 

The BRST Ward identity implies 13 



ik v (Av- ■ ■) 



r^ 4 + M£'/£ 



(124) 



where -A^ and <p are the truncated Green's function for 
emission of a gauge boson and Goldstone boson with mo- 
mentum k. This gives the equivalence theorem factor C 



C 



1 r, 



k 2 /i 



MTv> A + M£'/£ 



(125) 



Gauge bosons with polarization k 1 * in truncated Green's 
functions are replaced by —iMCip. Eq. (|123p gives C = 1 
at tree- level. 

The Ward identity Eq. (|124|) is exact. The scatter- 
ing rates are given by the S-matrix, which is given by 
multiplying the truncated Green's function evaluated at 
k 2 = A/p hys by the wavefunction renormalization factor 

V^a, and using the physical polarization M/E cor- 
rections are introduced by replacing the polarization cl 
by the momentum, = fc^/Mphys + Se 1 ^, and then ne- 
glecting the (5e^ term. The fc^/Mphys term can be com- 
puted using the Ward identity: 



(A L ...\S\---) 



-iC 



M 



-m 



1/2 



phys 



(£■ ■■)!*»= 



M- 



(126) 



The r.h.s. is evaluated at the physical gauge boson mass 
even in R^ gauge where the Goldstone bosons have mass 
M£7V£ We will write Eq. fTHfy as 



(A L --\S\-) 



•> 



k 2 =M 2 



C 



m m 1 / 2 

M phys ml/ 2 



(127) 



where fR^ is defined as though it is the wavefunction 
renormalization factor for tp, 



dr^ 



dfc 2 



(128) 



where the derivative is taken at a fixed (^-independent) 
reference momentum. The quantity we need, the 
l.h.s. of Eq. (|127p is independent of the choice of %K V . 



1/2 



fc 2 =A/ 2 



is a fake S-matrix element for 



ip, 



and defines the quantity on the r.h.s. of Eq. (|120p . It is 
computed like an S'-matrix element, but it is gauge de- 
pendent because it is not the scattering amplitude for a 
physical state. 14 The gauge dependence cancels that in 
£, so that the left-hand side is gauge independent. We 
have checked this by computing both quantities in 
gauge. The l.h.s. of Eq. (|127[) is also //-independent. On 
the r.h.s. of Eq. (|127p . £ and the (p fake S'-matrix element 
are both //-independent. 

The results tabulated in Sec. IXIVI give the S-matrix 
element for fermions and transverse gauge bosons. For 
longitudinal gauge bosons, we list the expression for the 

1/2 

Goldstone boson (p, which is why we have included 9^ 
so that it behaves as much as possible like an S-matrix 
element. The results in Table Hill are given in i?j=i gauge, 
and Vi v is evaluated at the reference momentum k 2 , = 



13 See Ref. |4l| for an excellent discussion of BRST Ward Identities 
and the Equivalence Theorem. 



14 This is true at one-loop even if /£, = 1 so that the Goldstone 
boson and gauge boson masses are equal at tree-level. 
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M 2 , which is the tree-level gauge and Goldstone boson 
mass. 

The one-loop results we will need are: 



M, 



phys 



M 



= 1 
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Is 



I log 



M 2 
M 2 



9 



" F 



2-z 
2 + z 



(129) 



where z = M^/M. This expression, and the ones below 
can be analytically continued to z > 2 using 

yfw tan -1 y/w <-> - 



ui tanh 




The expression for £ in 't Hooft-Feynman gauge in the 
MS scheme is 



47tt 73 1 

777 + 7; n F 



16tt 2 



4 12V3 12 3 



z 4 7z 2 5\, 



2 

y 



Uz 2 15 



4 y vi^ 2 



tan 



2-z 
2 + z 
(131) 



Notice the expression is /i independent as discussed ear- 
lier. 



XIV. RESULTS 

The results needed for a NLL analysis of hard scatter- 
ing processes in the SU{2) theory are summarized here. 
The computations for external (massless and massive) 
fermion and scalar states were already discussed in de- 
tail in Refs. @, H, 0]- They were presented previously 
using the sum-on-pairs form discussed in Sec. IX El They 
are tabulated here using the collinear+soft function form 
discussed in Scc. lXIl The new computations are those for 
external gauge bosons and Higgs bosons. The standard 
model results can be easily derived using the results pre- 
sented here; however, the expressions are considerably 
more involved than those for the SU (2) theory because 
custodial SU (2) is no longer a symmetry. Thus one has 
to include the effects of 7 — Z mixing, Mw ^ Mz, and 
the b — t mass difference. For this reason, the standard 
model results needed to compute all standard model hard 
scattering processes are given in a companion paper [46| . 




FIG. 26: The upper graph shows the topology of wavefunction 
corrections included in Table [TTT1 Topologies such as those in 
the lower graph are not included. 



The results for external transverse gauge bosons were 
computed by matching onto SCET operators where the 
gauge field is given by the SCET field, 



71, p 



iD% = id n + gA n , p , (132) 



The field does not produce longitudinal gauge bosons 
to leading order in the SCET power counting. The am- 
plitude for longitudinal gauge bosons is computed using 
the Goldstone boson equivalence theorem discussed in 
Sec.GTn] 



A. Collinear functions for the SU(2) theory 

The one-loop collinear functions are summarized in Ta- 
ble lllll The function fg is given in Appendix B of Ref. 0] • 
The tp row is for fermions, <fi for scalars, W± for external 
transversely polarized gauge bosons, H for the physical 
Higgs field, and ip a for the unphysical Goldstone bosons, 
which are used to compute longitudinally polarized gauge 
bosons using the equivalence theorem. The upper section 
of the table gives the results for fermions, scalars and 
bHQET fields, including wavefunction rcnormalization, 
assuming that the only wavefunction graphs are gauge 
boson corrections to the propagator with the topology 
shown in the upper graph in Fig. 1261 The lower section of 
the table gives the results for fermions, scalars, bHQET 
fields, transverse gauge bosons, the Higgs, and unphysi- 
cal Goldstone bosons, where the wavefunction contribu- 
tion has been separated out explicitly. 7^, etc. are the 
anomalous dimensions of the fields and i5£H is the finite 
wavefunction factor needed in the MS scheme, from the 
residue of the pole in the two-point Green's function func- 
tion, (1 + (5«K)/(p 2 -to 2 ). 



1. Running from fj,h to fii ~ Mz 

The collinear anomalous dimensions are given using 
the 7c column in Tabic HTT1 
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TABLE III: The collinear anomalous dimension and low-scale matching. Lm = log(Af 2 //i 2 ), L p = log(n ■ p)/fj,, and 7 = E/m. 
The rows are ip: fermion, <j> non-Higgs scalar multiplet, h v HQET field, B±: transverse gauge boson, h: Higgs, Lp a : Goldstone 
bosons (i.e. longitudinal gauge bosons using the equivalence theorem and multiplying by £). The results are in i?£=i gauge. 
yw,h,<p and 9iw,h,ip are the wavefunction contributions. 



For the fcrmions ip: 



^(4, g ^-3 



(133) 



with C F = 3/4. 

The gauge field anomalous dimension at one-loop in 
i?S=i gauge is 




FIG. 27: Collinear matching graphs for [VF^V]- The ® is the 
[W^] operator, the solid line is if) and the double line is . 
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(2C A - b ) 



(134) 



where bo is the coefficient of the first term in the /?- 
function, 



/' 



dg 
dfi 



-b 



16tH 



(135) 



so that the collinear factor 7 c for transverse gauge bosons 
is 

Wt (transverse W 1 * 2 ' 3 ): 



a 

47T 



(4CUL P - b ) 



(136) 



with bo = ll/3CA — 2/3np and Ca — 2. The final expres- 
sion is gauge invariant. There is a factor of the coupling 
constant g(nh) in the high-scale matching for each gauge 
boson. The — bo term in Eq. (|136[) converts that to g(m). 
The bo term is absent in the expression for f in Ref . [27| ■ 
The unphysical Goldstone bosons and Higgs are renor- 
malized together as elements of the scalar doublet (j>, 
<p,H: 



a n-p 
-C^4 log — 



(137) 



with C F = 3/4. 



2. Matching at m ~ Mz 

The matching corrections at /i/ ~ Mz have to be com- 
puted in the broken electroweak theory, using Table IIIII 
The collinear matching Dc depends on the mass spec- 
trum of the broken theory. The low-scale matching for 
the fcrmions is: 



exp D c Vjj 
expD^ E L 



(138) 



One can think of [W"r£] as a single collinear object which 
is matched in the effective theory onto the fields vl and 
El (see Fig. [27]). Custodial SU(2) invariance implies 
that vl and El get the same collinear matching. In 
the standard model, custodial SU(2) is broken, and the 
corresponding expression reads 



EW J 



(139) 



The collinear Wilson-line Wg W on the l.h.s. contains W 
and B fields, whereas on the r.h.s. Wl contains only pho- 
ton fields. The two Dc functions in Eq. Q139|) are no 
longer equal. One also has different collinear functions 
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Dq for the left and right-handed quarks, top quarks, etc. 
For this reason, the complete standard model expressions 
are deferred to Ref. [46j ] . 

The collinear matching for the scalar doublet is: 



exp D 



c 



^expZ^T^ H-j-expDi; 



(<t>-^<p) ,„3 

(140) 




where (p + , ip 3 have the same collinear matching by 
custodial-SE/(2) invariance. There is an interesting sub- 
tlety in the scalar doublet matching for the standard 
model, which is discussed in Ref. (46[. 

The matching for the transverse gauge bosons is: 



FIG. 28: Pair production q(pi)+q(j>2) — > q' (p4)+q' {p3). Time 
runs vertically. This also defines the variable for q + q — > g + g 
if particles 3, 4 are replaced by gauge bosons. 



[W*W£] 



(141) 



Again, custodial SU(2) ensures that all W a have the 
same collinear matching, and we do not have to worry 
about W 3 — B mixing. 
The Dc functions are: 
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where 



3q 

47T 

^L{2F w + 2f s (l,l)} + ±6K w 



(142) 



M w n-p 
2 log —f- log 



log 



1, 2M W 



M 2 



5^ 
~L2~ 



1 
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in terms of which, the soft anomalous dimension and low- 
scale matching arc 



is 



J(a(/i), U 



Tj U s {ni,nj) 



T 3 - ^(ni.tij-) 



(.145) 



For the Sudakov problem, there arc only two particles, 
and — m ■ 7i2 = 2, so that [7$ = 0, and the soft contribu- 
tion vanishes. The Sudakov computation is thus entirely 
given by the sum of the collinear functions for the two 
external legs. 

The soft function has a universal form when written 
in the operator form Eq. (|145[) . For numerical compu- 
tations, it is more convenient to choose a basis of gauge 
invariant operators, and write the soft-anomalous dimen- 
sion and matching as a matrix in the chosen basis. We 
compute the soft factor Y^Uj) Ti ■ TjUs(rii,rij) for some 
simple cases below for an SU(N) gauge theory, to give 
the well-known results for the soft anomalous dimension 
for QCD [86[. We need certain standard soft matrices 
which will then be used to compute the soft-running in 
the broken gauge theory. 



The wavefunction factors 9i V! H,W are those for the 
standard model with sin 2 8w — * 0, and can be found in 
Ref. [10, [H, Hl| . The A (tfcj)) 2 coupling enters the calcu- 
lation through the run dependence of the wavefunction 
factors. 



B. Universal Soft Function 



The universal soft functions is 



U s (ni,iij) = log 



-m ■ n.; — i 



0+ 



(144) 



Anomalous Dimension for qq 



qq 



As a sample computation, consider the soft function for 
qq — > qq for fermions in the fundamental representation 
of an SU (N) gauge theory. The kinematics for qq — ► q'q' 
is illustrated schematically in Fig. |2"81 where the incoming 
and outgoing particles have momenta p%, p2 and P3, pi, 
respectively, and we work in the limit s,t,n> M 2 3> m 2 . 

The 

2 



The external particles are all on-shell (pf 
Mandelstam variables are s = (p\ +P2) 2 , t 
and u = (p 3 - pi) 2 . 

An operator basis for gauge invariant operators is 



mf). 
(Pi ~ 



Pi 



Oi 



^T a i> 3 </> 2 T>i =T a ®T a 
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Invariant 


SU{2) 


SU(3) 


SU(N) 


Ca 


2 


3 


N 


Cf 

c d 

Ci 


3 
4 



3 
16 


4 
3 
5 
3 
2 
9 


N 2 -l 

2N 
N 2 -4 

N 
JV 2 -1 

4N' Z 



so that 



TABLE IV: Group theory invariants for SU(2), SU(3) and 
SU(N). 



Q 2 = tp4fp3 = 1 ® 1 



(146) 



where only the group theory structure has been shown. 
The group theory factors can be simplified using 

T a T a = C F 1 
T a T b T a = -h(Jj^j T b 

T a T b T a T b = C 1 l®l + ^{C d ~C A )T a ® T a 
T a T b T b T a = ^ 1 ^ 1 + I + C A ) T" ® T a 

(147) 

in the notation of Ref. [90]. The group invariants are 
listed in Table [IV] 

For scattering kinematics, s > 0, t < 0, and u < 0, and 
the variables T, U are defined by [86[ 

T = log h wr , 

s 



— u 

U = log h . 

s 



(148) 



The object that enters the soft running and matching 

is (Uij = U s (ni,nj)): 



(ij) 



iCdr-r + iCAra n 



Cm 







n 

where 



[Uu + U 23 ] - [U 13 + U 2i ] 

[Uia + U 23 ] + [U 13 + U 2i ] - 2 [U u + U 3i ] 



(149) 



Ui 2 


= U 34 = 


—in 


u 13 


= u 2i = 


log 

s 


U u 


= u 23 = 


log — 

s 




= 2(T- 


u) 


>2 


= 2(T + 


U) 



(150) 



-2i7rCVl 



±C d (T-£/) + ±C A (T + [/) 2(T-U) 
2d(T-U) 



(151) 



which agrees with VA and VC of Ref. 0, and Ref. gi] 

We will need the soft matrix Eq. (|151[) when the fields 
arc in the fundamental representation of SU(3), 



6 0) = _^ Ti . Tj . u 8 (ni,i 

(ij) 



--wrl 



|T+|C/ 2{T-U) 
4 







(152) 



and when they are in the fundamental representation of 
SU(2), 

6 (2) = -J^ti-tj C/sK,nj) 

(ij) 



-int 



(T + U) 2{T-U) 
HT-U) 



(153) 



If two of the fermions (say 1, 2) are in the fundamental 
of SU(N), and the other two are singlets, then the only 
allowed operator has the form 1 (g> 1, and the soft factor 
is 



J^Ti-Tj U s (ni,rij) = -iirC F 

(ij) 



with values 



6 (3)/ = _^ 

3 



6< 2 >' = --i7r, 
4 



(154) 



(155) 



for SU(3) and 5(7(2), respectively. If all four fields are 
singlets, the soft factor vanishes. 
For a U(l) gauge theory: 

© (1) (<?1,92,93,94) = -^Tj ■ Tj U s (n, t ,nj) 

m 

= 9l 92^12 + 9394 ^34 + 9194^14 

+9293 ^23 ~ 9193^13 ~ 9294^24 
• 7r / 2 , 2 1 2 , 2\ 

= (9i +9 2 + 93 + 9 4 ) 

+ (9i94 + 929 3 ) T - (9193 + 9294) J/ 

(156) 

which agrees with VIIB of Ref. and VIIC. In most 
cases, we need this result for q\ — q 2 = qi and q 3 = 54 = 
9/, 

s(1) (9/>9*) = -z2<li- Ij U s {ni,rij) 
(ij) 



■iir (9- + 9/) + 29,9/ (T - V) 



(157) 
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FIG. 29: Quark scattering q(pi) + q {ps) — > qipz) + q'ipi)- 
Time runs vertically. This also defines the variable for q+g — > 
q + g if particles 3, 4 are replaced by gauge bosons. 



2. Anomalous Dimension for qq qq 

The scattering kinematics is given in Fig. [29l and the 
kinematic variables are defined by s = (pi +P2) 2 , t = 
(pi — Pi) 2 and u = (p^ — Pi) 2 - For scattering, qq — ► qq, 
the soft matrix in the i-channel basis Eq. (| 146)) becomes 

-Y,Ti-Tj Us(ni, nj ) 

" \{C d + C A )U -C A T 2U 
2dU 

(158) 



2C F (T + 



which agrees with VB of Ref . Q ■ 



3. Anomalous Dimension for Heavy quark pair production 

The soft matrix is identical to that for qq — > qq, and 
agrees with VC of Ref. 0. 



4- Anomalous Dimension for qq — > gg 

The scattering kinematics is given by Fig. [55] with s 
(Pi + P2) 2 , t = (p 4 -pi) 2 and u = (p 3 - pi) 2 
The operator basis for — > is 



Or = 1 ® <5 AB = miAfAf 

2 = T c ® d ABC = d ABC q 2 T c qi A A A B 

3 = T c ® z/^ sc * = if ABG q2T c q\A±A B 

The soft operator in this basis is 

(ij) 

= -in(C F + C A )l 





(159) 



U -T 



iC A (T + C/) §CU(tf-T) 
2(U-T) \C d (U-T) \C A {T + U) 



(160) 



which is the result of Ref. [86j . 

For QCD, this gives the soft matrix 



6 



(3,9) 



13 ■ „ 
vk\ 

3 



U-T 

U T + U ) U U ~ T ) 

2(U-T) UU -T) UT + U) 



(161) 



For SU(2), the d-symbol vanishes, so the operator ba- 
sis reduces to 



Oi 
2 



9291 WW 



(162) 



and the soft matrix is given by using the SU (2) values for 
the group constants, and dropping the second row and 
column, 



g(2,9) 



11 



-inl 



U-T 
2{U-T) (T + U) 



.(163) 



5. Matching for Fermion Scattering 

The soft contribution to the low-scale matching can 
be computed using the above results. Consider, first the 
case of a four-quark operator, where all the fields are 
SU(2) doublets, for example 



(164) 



where the index is 1 for t a ® t a and 2 for 1 <E> 1 in SU(2), 
The group theory sum needed for the soft anomalous 
dimension matrix is 

© 2 = -J^ti-t,- U s (n i ,n j ) = eW (165) 
(H) 

where &^ is the standard form Eq. (| 1 53[) . The soft 
anomalous dimension is 



IS = -©2' 
7T 



(166) 



At the low scale fii ~ mz, the operators Eq. (| 164|) match 
onto a linear combination of 

O22 = [^ L4 l a ^ L3 \[eL2l a e L i} 

O32 = [P- / L4JaHL3}[VeL2j°'VeLl} 

O42 = [Mi47aMi3][eL27 Q e^ 1 ] 

O52 = [P-L4jaV t i L3 ][i'eL2j a eLl) 



a 



(.2 



li4 7aML3][eL27 Qjy eLl] 



(167) 
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with coefficients C,-,- . The matching matrix is 



(168) 



At tree-level R is 



(169) 



The one-loop soft matching due to W exchange is given 
using Eq. flST)) 



R 



(i) 
s,w 



— 2\og—f- 

47T /i Z 



Eft 



Using Eq. (|165[) . the one- loop soft matching is 



R 



(i) 

sr,w 



^21og^ 

47T 



(170) 



(171) 



1? 



(i) 

Sf.W 



gam, one 



and the total soft matching is i? 1 - ** 
can see the great simplification due to custodial SU{2) 
invariance. Since W 1 ' 2 ' 3 all have the same mass, the soft 
matching depends on the ST/ (2) invariant combination 
t j • tj , so that the soft matching can be written in terms 
of the same matrix 6 2 that enters the soft anomalous di- 
mension. In the standard model, the sum has to be bro- 
ken into the charged W and Z contributions, and takes 
a more complicated form given in Ref. [46j |. 



6. Matching for Electroweak Gauge Boson Pair Production 

The kinematics for the electroweak gauge boson pair- 
production process is shown in Fig. 1301 The operator 
basis is 



2 



3 ! 



(172) 



where only the gauge structure has been shown. Note 
that e ahc W^W\ ^ since the two W fields have momen- 
tum labels P3 and P4 which arc different. Each gauge 
structure has 10 possible momentum structures, which 
are given in Ref. [3 11 ] - The EFT running and matching 
is the same for all 10 amplitudes, so the momentum de- 
pendence can be ignored for our analysis. 
In the basis Eq. (fT72"l) 



e 2 = -J> 



tj Us{n i ,n j ) = e^ (173) 



w(po 




W + (p 3 ) 



FIG. 30: Pair production 5 (pi) + 5(252) 
Time runs vertically. 



W + { P3 ) + W-{pi). 



where <S^ 2 ^ is given in Eq. (|163[) . The soft anomalous 
dimension between the scales Q and m ~ Mz is given by 
Eq. (fT66|) with Eq. (fl73j) for 6 2 . 

At the low-scale /i; ~ Mz, the operators Eq. (1 1721) 
match onto 





= ^L2^LlW^W 3 


60 


= V L2 V LX W^W^ 


63 


= VL2VLlW\Wl 


64 


= e L2 e L1 W+Wi 


Ch 


= e L2 e L iW^W+ 


O e 


= e L2 e L1 WiWi 


O7 


= y^LiW^Wl 


8 


= VL2ZLlWlW+ 


<5 9 


= emvhiWlW^ 


610 





(174) 

where W ± = {W 1 + iW 2 )/\/2. The subscripts 3,4 rep- 
resent outgoing label momenta p 3 and 754, and the gauge 
indices are to be treated as those on a quantum field, 
i.e. they represent the charge on the annihilation opera- 
tor. These operators are used like terms in a Lagrangian. 
Thus ee — ► W + (ki)W~ (k 2 ) is given by C4 with p 4 = k\ 
and P3 = k 2 , plus C5 with p^ = ki and p$ = k 2 . 
The tree-level matching is 



Ci = 





1 

1 

71 



(175) 
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The one-loop soft matching due to W exchange is given 
by Eq. (flTT) using the value of 6 2 in Eq. (fT73"|> . 



The tree-level matching is 



7. Matching for qq — > (pip 



For longitudinal W production, we also need the re- 
sults for external unphysical Goldstone boson tp holds, 
which arc contained in the Higgs multiplet cj>. The oper- 
ators are 



(h 



(176) 



The gauge current i T a D ^ — $ T a (j>) produces 
operators of this form, weighted by a label momentum 
V(f, which is included in the operator coeffi- 



factor 7^4 



cients, and is antisymmetric in 3^4. 

The group theory sums needed for the soft anomalous 
dimension matrix are 



© 2 = - *« • *j Us ^ n *) = s(2) • ( 17? ) 

(ij) 
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2V2 



(179) 



and the one-loop matching is given by Eq. (|17ip using 
the values of and & 2 for the current case. 



C. Final Result 

The final result for the scattering amplitude is then 
given by using Eq. (|T4[) . with the running and matching 
given by the sum of the collinear and soft contributions, 
as in Eq. (fT07|l . 



D. Comparison with Previous Results 



At the low scale, the operators match onto 



Ox 
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= v L2 vliHhp\ 


64 


= v L2 vlw\H z 


65 


= VL 2 VLiHiH 3 
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= e L2 e L1 ip^ip^ 


67 


= ez^eLivtVi 


Os 


= e L2 eLii?4yi 
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= e L2 e L iiplH 3 


610 


= ZLlZL\ti±H?, 


On 


= v L2 eLif\ft 


612 


= v L2 e L iH 4 ip+ 


613 


= e L2 u L iip^ipl 


Ol4 


= e L2 v L iip^H 3 



The clectroweak radiative corrections using SCET were 
computed previously in Refs. 0, 0, Q for processes not 
involving external gauge bosons. The results given here 
are much more compact, the hnal results are the collinear 
functions for each particle given in Sec. IXIV Al and the 
soft functions computed in Sec. IXIV Bl The previous re- 
sults were given in terms of kinematic variables such as 
Q 2 or s, t and u. and summing over pairs of particles. 
The collinear and soft functions given here are more con- 
veniently written in terms of hi ■ pi , and rii ■ rij , which 
depends on the choice of n. The n dependence cancels in 
the hnal result between the collinear and soft contribu- 
tions. The kinematic conversions needed to compare the 
two forms of the result are given below. 

In the Sudakov problem, Q 2 = — (pi +p 2 ) 2 , where the 
two momenta are incoming. One can rewrite 



log— = log 



ni-pi n 2 -p 2 

1- log (180) 



(178) 



where —hi ■ h 2 = 2, since n\ = (1, n) and rio = — (1, — n). 
Using this, the results in Table I of Ref. Q and Table I 
of Ref. Q can be obtained by adding pairs of rows in 
Table HU 
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Tabic III of Rcf. [6[ for the matching at heavy particle 
thresholds is equivalent to Table |TTJ If particle 1 is light 
and particle 2 is heavy, the relation 



O 2 1 m\ 
l0 V~2 1Og ^ 



log h log ■ 



= log 



fl\ - Pi 



m 2 

log(2 72 ) (181) 



using fi2 ■ V2 = 272 can be used to show the equivalence 
of the two forms. 

Tabic IV of Ref. [|| for the heavy-heavy case is also 
obtained from Table ITTl using 

log(2u>) = log(2 7 i) + log(2 72 ) (182) 

for w = Vi ■ V2 3> 1. 

XV. SUDAKOV FORM FACTOR 

The SCET results are applied to the Sudakov form fac- 
tor. This is an instructive example, and allows us to con- 
nect the SCET formulation with the results of Magnea 
and Sterman [Z^] (and more recent work in Refs. [9ll.l92l|) 
obtained using factorization methods. The SCET de- 
composition of the amplitude into collinear and soft is 
related to, but not exactly the same as, the factorization 
decomposition, and it is useful to study the difference. 
In addition, we can compare the Sudakov form factor in 
a theory with massive gauge bosons with the Magnea- 
Sterman analysis for QCD. 

The spacclikc Sudakov form factor F(Q 2 ) is defined as 
the particle scattering amplitude by an external current, 
with momentum transfer Q 2 = —q 2 > 0. It is conve- 
nient to compute the form factor in the Breit frame (see 
Fig. [T5)) . where the particle is back-scattered, and the 
momentum transfer q has q° = 0. One can compute the 
Sudakov form factor in this frame using Eq. (|5§|) with 
r = 2, where particle 1 and 2 have the same quantum 
numbers, and n\ = —fi2, ri2 = —n\. There is only one 
gauge invariant amplitude, so there is no matrix struc- 
ture to the amplitude. We work in the large momen- 
tum transfer limit, so that M 2 /Q 2 power corrections are 
dropped. 

The form factors for fermion scattering by the vector 
current j M = ^7^7, the scalar current js = ipip and the 
tensor current jl^f = ipcr**" ij) are defined as Fi/(Q 2 , M, /x), 
Fs(Q 2 ,M, fi) and F T (Q 2 , M, fx), where the operators are 
renormalized in the full theory at a scale fi. 15 They are 
normalized to one at tree-level. The anomalous dimen- 
sions of the operators in the full theory are 



-lv{a)f = 0, 



15 We also give results for gauge scattering by the field strength 
tensor F 2 = F^F^" in Appendix fA] This provides a check on 
the collinear anomalous dimension for the gauge field. 



d/x 

u—tpa^ib 
d/z 



-7s(a) ipil) = 6—C F ipip , 

- ot — 

-7t(«) tpa^ip = -2—C F ipa^i), 
in 



(183) 



where we have given the one- loop values, and defined our 
sign convention for 7. The form factors must satisfy the 
equations 



p,—F {Q 2 ,M^) 
d/j, 



-7o(aM) F Q (Q 2 ,M,n). 



(184) 



In SCET, the form factor is given by running the oper- 
ator from fi to fih using the full theory anomalous dimen- 
sion 7/, matching onto SCET at p,h, running in SCET 
from fjLh to fii, and then evaluating the matrix element 
at fii. The result is more conveniently written for logFb, 
where = V,S, T: 

1 du 2 
\nF (Q 2 ,M,^) = - ^70 (a^)) 

+C (a (n h ) , log % 



^ du 2 r o 2 

r(a(/i))ln^+B(a(A)) 



f> 2 



M 2 \ ( M 2 \ Q 2 

-D ( a (m) , log — + L>i la (m) , log — In — 
M; / V M; / ^1 

(185) 



where 



b = vti + n 2 



(186) 



from Eq. (|68p . Explicit one-loop expressions for the form- 
factor and the different terms in Eq. (| 185[) can be found 
in Appendix [Al 

Note that there is a single-log term in the low-scale 
matrix element. The EFT quantities T, B, D0.1 do not 
depend on the choice of fermion current, but the full 
theory anomalous dimension 70 and high-scale matching 
Co do. Eq. (|185|) is independent of [ih and fj-i, which 
implies the consistency conditions 



lo{a(p)) 



d ( / \ M 2 



d/z 



2D X a(/i),log 



-H—Dq a(/z),log — 
dfi \ \i A 

-r(a(/i))ln^ +J B(a(/i)) 

-A*x~C(a(M)><3/M) ■ 
d/j, 



(187) 
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Eqs. (|185I187[) arc valid to all orders in perturbation the- 
ory. 

Differentiating Fq(Q 2 , M, //) w.r.t. logQ 2 gives 



d In Fp{Q 2 ,M, fi) 
SlogQ 2 



M 2 

K ( a(/i),log-^- 



+G a(/i),log- 



/'- 



,(188) 



where -Ko is the logM 2 / (x 2 dependent part of the deriva- 
tive, and Go is the rest. The form Eq. (|188p is highly 
non-trivial, because the logM 2 and logQ 2 dependence 
has been separated into different terms, as shown below. 
Using Eq. (|185[) . one gets 



(2) 

D\ for an SU{2) gauge theory with degenerate Higgs 
and gauge bosons is @ from the results of Ref. [l5| 



Df\jjL = M) = C F n F T F 



-5C(3) 



47T 2 112 

~ + ~27 
15V3?r 391 
4 18" 
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V3Cl 2 (7r/3) 



(192) 



D\ depends on the full mass spectrum of the theory. 
Differentiating the form-factor evolution equation 



d/z 



-7o(a(/z)) 



(193) 



w.r.t. logQ 2 gives [42 



M 2 

G = 2Di ( a(/i(),log— 2- 



^a(Mh),log- 



<91nQ 2 

^ d£t 2 



(189) 



The r.h.s. of Eq. (|189[) is independent of and hi due 
to the consistency conditions Eq. (|187[) , and can be writ- 
ten in terms of a(n) and Lm = log M 2 /fi 2 using the (3- 
function to transform all coupling constants to the scale 
/!. Equivalcntly, one can choose \ih = jjLi = ji in Eq. (|189|l . 



M 2 

R'o + Go = 2D, [ a ( M ),log — 



flln Co (oO*), log ^) 



so that 



K 



Go 



ainQ 2 



M 2 

2D i ( a(/i),log — 



d In Co (a(/i),logHr 



, (190) 



<91nQ 2 



(191) 



Z?i and hence if are independent of the operator O, so 
we can drop the subscript on K . Eq. (|19ip justifies the 
form used in Eq. (|188p . Eq. (|191[) follows from the low 
scale matching which depends on M. D + D\ log Q 2 / [if 
being linear in logQ 2 , as proved in Rcfs. [|36|, l66f , The 
derivative of the low-scale matching w.r.t. logQ 2 is then 
independent of Q 2 , and gives the expression for K. 

The low scale single-log matching D\ is independent 
of the operator O, and is universal. D\ is the universal 
function of Ravindran, Smith and van Neerven (93l . 
Q for a massive gauge theory. D\ is related to the cusp 
anomalous dimension by Eq. (|187[) . and so should have 
the maximally non-Abclian color structure of the cusp 
anomalous dimension. This is true at two-loops, where 



(194) 



This is clearly satisfied by Eq. (|189j) since fi does not ap- 
pear on the r.h.s. However, K and Go separately depend 
on [i when the r.h.s. is written in terms of a(p), as in 
Eq. (fT9T|) . Eq. fl94|l gives 



d/Lt 

V-r-Go 
d/z 



-j K (a (/i)) 
Ik {a ill)) , 



(195) 



since a(fi) is the only common variable to both K and 
Go- Jk does not depend on O or Q 2 , since K does not, 
and it does not depend on M 2 since Gn does not. Using 
K in Eq. (MD and Eq. flgg]) gives Hf 



7 x(a(/x)) = 2I>( M )) 



(196) 



which shows that the single-log term in the SCET anoma- 
lous dimension is related to the cusp anomalous di- 
mension. Eq. (|195|) also follows from Eq. (|187[) using 
Eq. (TITO . 

Integrating Eq. (|188p gives 



In F G (Q 2 , M, //) = \nF {Ql,M,fi) 

1 z^ 2 d^ 

2 4 £ 2 



M 2 

K ( a(/x),log — 



e 



Go ( a(M) 5 l°g ^ 2 



(197) 



and integrating Eq. (|195[> gives 
G o (a(0,0) = Go Ca(M), log ^ 



1 /* d M 2 

(198) 



so that 



lnF (Q 2 ,M,fi) = lnF (Q 2 ,M, M ) 
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e 



'Ql 

+G (a(O,0) 



M 2 

K ( a(»,log-^- 



o / — 7if(«W) 



(199) 



which is the expression for the Sudakov form factor in a 
massive gauge theory written as in Ref. [42| . This equa- 
tion is equivalent to the SCET form Eq. (|185l) ; a detailed 
comparison can be found in Sec. IV of Ref. [y]. 



A. Massless QCD 

The entire analysis above can be repeated for a mass- 
less gauge theory such as QCD. In Eq. (|185[) . the high 
scale running, high scale matching and running between 
Hh and [ii using the SCET anomalous dimension are un- 
changed by the infrared structure of the theory, and are 
identical in the broken and unbroken gauge theory. The 
only difference is the low scale matrix element D. In 
QCD, the low-scale matrix element D is infrared diver- 
gent, and has the form 

D (a(jn),e) + Di(a(jit),e)\n% (200) 

where the 1/e terms are infrared divergent, and are not 
cancelled by renormalization counterterms. So all the 
results of the previous section hold with the replacement 
of a functional dependence on log M 2 / fx 2 by one on 1/ e. 16 
Explicit one-loop expressions for the massless form factor 
are given in Appendix [X] 

The form-factor computed for a massless gauge theory 
in pure dimensional regularization is IR divergent, and 
has the functional form 



(201) 



F ( a s (fi),e, -jj 



The explicit fj, dependence in Fy(Q 2 , /J, 2 ) (i.e. for the vec- 
tor form factor) is cancelled by the implicit [i dependence 
in a, so that Fy is /i-indcpendcnt. The EFT computa- 
tion splits the form factor into a high scale part which 
depends on Q, but is IR finite, and low-scale matrix ele- 
ments -Do,i which are IR divergent, but Q independent. 
The total low-scale matrix element has the linear-in-log 
form Eq. 1200]) . 

The Magnea-Sterman expression for the Sudakov form 
factor in massless QCD is given by Eq. (|199[) with the re- 
placement K — > K(a([i, e)). Magnea and Stcrman stud- 
ied only the vector form-factor, for which Fv(Qq = 0) = 



16 This does not mean that the functions are identical. Using M 
as an IR regulator instead of e can change the subleading diver- 
gences and finite parts. Compare Eq. (I A6t and J A 10 I t 



1, since it is the forward matrix element of the quark 
number current. They dropped In F(Qq) in Eq. (|199[) and 
set Qo — > in the limits of integration. We have used 
SCET in our computation, and been able to compute 
the form factor for any current. We cannot take the limit 
Qo — > 0, since by construction, the SCET results are only 
valid for large Q, and power corrections are neglected. 
In the massive gauge theory, this is because we have ex- 
plicitly dropped M 2 /Q 2 power corrections. In massless 
QCD in perturbation theory, it might seem that there is 
no other scale in the problem, so we could use SCET for 
Q — > 0. However, this apparent scale independence is il- 
lusory; perturbation theory knows about the scale Aqcd 
through the coupling constant, and we have neglected 
non-perturbativc Aq CD /Q 2 corrections from higher di- 
mension operators. 



XVI. CONCLUSION 

We have extended our previous results on high-energy 
scattering to include external longitudinally and trans- 
versely polarized Higgs bosons. This allows us to com- 
pute any high energy scattering process in the standard 
model. 

The factorization structure of the effective theory al- 
lows us to write the results in terms of single-particle 
collincar functions which only depend on the energy of 
each particle and a universal soft function that only de- 
pends on the particle directions. This holds for both the 
anomalous dimension and the finite matrix clement in 
the broken gauge theory. 

The electroweak corrections are substantial and in- 
crease rapidly with energy. They are largest for pro- 
cesses involving transversely polarized electroweak gauge 
bosons. The one-loop electroweak corrections are more 
important than the two-loop QCD corrections, in the 
cases we have studied. The EFT method allows us to 
compute the hard cross-section with a theoretical uncer- 
tainty of better than 1%, except for transverse W pro- 
duction, where the uncertainty is < 10%. 

The EFT result can be used to give the renormalization 
group improved cross-section if the high-scale matching 
is known. The high-scale matching is two powers of L 
smaller than the other corrections in the log-counting. 
Numerically, the largest high-scale matching we have 
found (including QCD) is about 10% (for Wt produc- 
tion), but this is also the case where the other EFT cor- 
rections are large, and the purely electroweak corrections 
are several times the matching even at 2 TeV. As the en- 
ergy increases, the matching corrections decrease slowly, 
whereas the EFT corrections increase very rapidly. One 
can compute the cross-section in cases where the high- 
scale matching corrections are not known by simply us- 
ing the tree-level values. This introduces an uncertainty 
< 10% in the rate, but now we can compute the radiative 
corrections to any standard model process with an arbi- 
trary number of external particles. The collincar correc- 
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tions are given by simply multiplying the contributions 
from each external particle. The only tedious computa- 
tion is the soft group theory sum over pairs of particles, 
but this is easily automated. 

The EFT results neglect M 2 /Q 2 power corrections. 
The largest power corrections are M§/Q 2 power correc- 
tions to the electroweak contributions. The Aq CD /(3 2 
corrections to to the QCD contribution are much smaller. 
As discussed in Sec. [V] the unknown power corrections 
are only in the one-loop terms; the tree-level power cor- 
rections can be included using the exact tree-level am- 
plitude and phase space. Thus the EFT results can be 
used close to threshold. For the Sudakov form factor, the 
power corrections are 2% exactly at threshold. 

The collinear functions for all standard model parti- 
cles, and the soft functions for some important 2^2 
scattering processes, are given in Ref. [46[. The results 
in this paper sum the NLL series and the largest NNLL 
terms. The missing terms for a complete NNLL compu- 
tation are the scalar pieces in the three-loop cusp and 
two-loop non-cusp anomalous dimension, which can be 
computed in an unbroken gauge theory, and the two-loop 
\ogQ 2 /fi 2 term Dl in the low-scale matching, which de- 
pends on the mass spectrum of the theory. The two-loop 
matching for Wt production is required to reduce the fih 
uncertainty in the cross-section. 

The partonic cross-sections computed here have to be 
convoluted at a factorization scale fif with parton dis- 
tribution functions and possibly jet and soft functions to 
get hadronic cross-sections. This is a non-trivial task, 
but has been well-studied and is independent of the com- 
putations given in this paper. 



APPENDIX A: THE SUDAKOV FORM FACTOR 

We summarize the one-loop expressions for O = 
V, S, T, and the high scale matching for O = F 2 here. 

1. Massive Gauge Theory 



The form factors for O — V.S.T are 



F v 
F s 
Ft 



aC F 

47T 

a Cf 

47T 

a Cf 

47T 



-Q/M 



3L 



-Q/M 



-Q/M 



2 -3L 



7 2tt 2 

Q/M - g ^~ 

5 2^ 

7 2tt 2 



m t- - 



3Lq/m + Lq - - - 



(Al) 



The full theory anomalous dimensions are: 

tv(«(m)) = o, 



Is {a{p)) 
IT {a {n)) 



-6- 



4tt 



-Cf , 



2 —° F - 



(A2) 



The high-scale matching is: 



CM = ^c F 
c sM - °-£c F 



Lq 2 + 3Lq 



-Lq' 



-q 



7T 

Y- 2 



4L Q + ^ 



(A3) 



The SCET anomalous dimension is: 



B(a(ji)) 



C F = -A (a (/z)) 



47T 

-6^C F 

4-7T 



(A4) 



The low-scale matrix element is: 

/ M 2 
D I a (m) , log — T 



M' 



D l [ a ( W ), log I In- , 
M; / H>i 



Q' 2 



(A5) 



with 



D 



D! = 



a(m) C F 
Ait 

a(m)C F 

47T 



M 



r2 



AI 



r2 



- log — - 3 log - , 

fif fif 2 

21og — 



5tt 2 



(A6) 



The effective theory form factors are all equal to 



F EFT (Q 2 ,M,/z) 



i(fi) C f 



-in 



M 2 Q 2 M 2 
-2 log — log — - 3 log — , 

fj, z jir 2 



log — 

[i 2 

9 5tt 2 

" ~6~ 



(A7) 



and differ from the full theory form-factors by the high- 
scale matching Eq. (|A3[) . 



2. Massless Gauge Theory 

The form factors are 



F v (Q 2 ,fi) = 1 



4-7T 



C f 



2 2,3 



7T 



2 



L q z + 3Lq + — -8 
6 



2 2, 

-3 + T L Q 



2 2, 

-3 + T L Q 



-Lo 2 +4L + ^--8 



(A8) 
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Note that the 1/e and Lq terms do not have the same co- 
efficient for Fg and Fx since there is also a 1/e UV diver- 
gence that is cancelled by the renormalization countert- 
erm, and leads to the anomalous dimensions in Eq. (|183D . 

The full theory anomalous dimension, high-scale 
matching, and SCET anomalous dimension have the 
same values as for the massive case, Eqs. (|A2IA3[ IA4I) . 

The low-scale matrix clement is: 



Q2 

D (a (fii) ,e)+Di (a (m) , e) In — , 



(A9) 



with 



47T 



D l = —C F 
4ir 



2 

72 



(A10) 



Z?o.i can be computed directly using the EFT. The 
EFT graphs are scaleless, and vanish, so Z?o.i are given 
purely by the counterterm contributions which cancel the 
UV divergences in the effective theory, so that Dqi con- 
tain purely IR singularities in 1/e. 

The effective theory form-factors arc all equal to 



F(Q 2 ,e, l x) 



«(m) 

47T 



2 2 Q 2 3 
— + - log — - - 
e z e fi z e 



(All) 



and differ from the full theory form-factors by the high- 
scale matching Eq. (|A3[) . 



F 1 



For gauge boson scattering by the field-strength tensor, 
the anomalous dimension in the full theory is 



^F 2 = - lF 2(a)F 2 , 

90(9) P(9) « 9 , 

7f 2 = —5 = -7- 2o 

dg g 4?r 



and the high-scale matching is 



C F 2 (/z) 



4tt 



C A 



(A12) 



(A13) 



The last result gives a consistency check on the — bo term 
in the collinear anomalous dimension for gauge fields, 
using Eq. (fl87|) . 



APPENDIX B: COLOR IDENTITIES 



The operator is a color singlet, so 
£Tf = 0. 



(Bl) 



Then 



= (X>) - (p T ,) 



J2 m ■ t* + E (fi + u ) T ■ T i > (B2) 



so that 



^(/ i + /,-)T i -T j = -^/iT^Ti. (B3) 



<</> 



Similarly, 




(B4) 



so that 



V] Tj • Tj = -Ti • T.; (no sum on i) . (B5) 



APPENDIX C: SYMMETRIC FORM FOR THE 
AMPLITUDE 

We follow the argument below Eq. (f52")l . but starting 
with the initial form Eq. (|79p . In Eq. (fT9"|) , treat Jj as a 
function of Z(i) = log(A.;/^ 2 ) = log(nj ■p i )/n + logS i /^i, 
and S^rij • rij) as a function of Xij = log(fii • n,) — 
log((5j//Lt) — log(5j//i). Then Si independence of the total 
amplitude gives 

T,T,«^T,. Tj «. (CI) 



dZ, 



dXi. 



Differentiating w.r.t. 6j, j =/= i gives 

d Sij(Xij) 







dx 2 ' 



(C2) 



so that S is linear in X^ , 

= J(a00, Af, //)^ - g'(a(fi),M, fj,) . (C3) 
Eq. (jCip then implies 
51,(2,) 



T, : ■ T, 



^ T, ■ TjJ(a(fx), M, n) 



= -T,-T 4 J(a(/i),M,/x), (C4) 

using the relations in Sec. [B] so that 

Ii(a(fi), M, Ai,mj,/x) = —J(a(fj,),M, (i)Z(i) 

+E(a(/i),M, mi,fx) 

(C5) 
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and the total amplitude is 



log 971 = J(a(fi),M,fi) 



J2 T * ■ T * lo s 



n% - Pi 



^T.-T.logK-n,; 



Ti ' Ti 



E{a(fx), M, m t ,fi) + -g(a{fx), M, fx) 



(C6) 

using Eq. (|B3|) . It is convenient to replace n.^ ■ rij — > 
(— rii ■ rij)/2, and absorb the log(— 2) term in g so that 



logon = J(a(fx),M,fx) 



]T T, • T, log 



rii - pi 



/' 



T,-T,\n A —i 



(ij) 



E(a(fx),M, mi,(i) + -g{a{fx),M, fx) 



(C7) 



where the i0 + is given by analytically continuing ■ rij 
for spacelike momentum transfer to timclike momentum 
transfer. 

This is the form of the general amplitude, Eq. ([55]). and 
as promised, is written in a manifestly symmetric form 
in the particle labels. The only term not proportional to 
the unit matrix is 



J(«0), M, fx) Ti-Tj log 

(ij) 



-Hi ■ Uj 



-«0+ 



(C8) 



For the other terms, one can use • T 4 = Cptl to write 
them in terms of the color Casimir. 
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